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ABSTRACT 

The design of airplanes, ships, automobiles, and so-called ''sculptured 
parts" involves the design, delineation, and mathematical description of 
bounding surfaces. A method is described which makes possible the description 
of free-form doubly curved surfaces of a very general kind. An extension of 
these ideas to hyper-surfaces in higher dimensional spaces is also indicated. 

This surface technique has been specifically devised for use in the 
Computer-Aided Design Project at M.I.T. , and has already been successfully 
implemented here and elsewhere. 
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SECTION I 



INTRODUCTION 



The purpose of this work is to present the mathematics of a certain class 
of surfaces which are suitable for the design and description of arbitrary shapes. 
In the past, the subject of surface mathematics has been investigated, in ana- 
lytical geometry and in differential geometry, from the standpoint of the analysis 
of geometric properties of surfaces that already exist, but very little literature 
has been produced on the subject of the creation of such surfaces. As a typical 
example, the design of the hull of a racing yacht requires the description of a 
surface of considerable subtlety and complexity, and the process is traditionally 
carried out by purely graphical procedures which are exceedingly laborious, 
since they entail a large amount of trial and error Iteration in order to assure 
that the surface is completely described, and is smooth and "fair. " The design 
of automobile bodies and airplane fuselages is similarly tedious and time con- 
suming, although mathematical techniques have been applied to aircraft design 
for a number of years. 

A few papers have been written on the subject of fitting existing ship hull 
shapes by means of various types of polynominals, with the two-fold purpose 
of smoothing and interpolating the information contained in preliminary graphi- 
cally derived hull lines, and of replacing this graphical information with for- 
mulas and equations that will permit further analytical techniques to be applied, 
such as structural analysis and the discipline of fluid mechanics. But these 
mathematical techniques are applicable only when the surface has already been 
designed to some degree of completeness, so as to furnish enough information 
to make the mathematics work. 

The mathematical structure of the surfaces to be described in the follow- 
ing discussion has been devised to implement the surface design process itself, 
so as to make it, from the designer's standpoint, extremely natural and easy. 
The designer himself need not know or care about these internal mathematical 
details, any more than he needs to know the specific composition of the pencils 
with which he writes or the mechanics of the splines with which he now draws 
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curves. The mathematics is relatively simple, but it is nevertheless too com- 
plicated for hand calculation, and is designed for use on a computer. 

In the design of a three-dimensional object , whether it bean airplane fuse- 
lage, an automobile body, a ship's hull, or a single sculptured part of a machine, 
toe designer requires a system which will permit him to define a surface with a 
minimum of input information, and then to modify this surface, if he feels so 
inclined, either by changing toe original input, or by adding more design con- 
straints to the system. 

As a specific example, suppose a designer wishes to design an airplane 
fuselage, using the SKETCHPAD system. ^' ^' ^ He would like to be able to draw 
the outline of flie airplane as seen from the side, the ouUine of the airplane as 
seen from above, and some arbitrarily selected section midships. With these 
three arbitrary curves designed, he would like to have the computer automati- 
cally and immediately generate a "fair" surface and display this surface to him 
m sufficient detail so that he could m^e appropriate judgments. If Ae surface 
so generated does not satisfy him, he would perhaps like to modify his original 
design curves, or else he might perhaps like to add other new sections and have 
the computer automatically and Instantiy re-fair the surface to fit this additional 
information. 

The following sections describe a very simple, flexible and general class 
of surfaces which are able to fulfill these requirements. It will be shown that 
a single algorithmic structure and essentially only two symbol types serve to 
provide the following features: 

1. Smooth, fair surfaces can be defined by a minimum number of curves, 
and then adjacent surfaces can be designed to match position, slope, 
curvature, and indeed any desired order of derivative along the ad- 
joining boundaries. 

2. The design curves that define the surface can be of any kind whatso- 
ever, including circles, second-degree curves, polynominals , trans- 
cendentals, and also sketched curves with no known mathematical 
formula whatsoever. 

3. Some classic surfaces are not necessarily members of the family 
of surfaces to be described; nevertheless, these classic surfaces 
can be matched along their boundaries to any order of derivative 
desired. 
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INTRODUCTION 3 

4. The arithmetic involved in constructing these surfaces is extremely 
simple and, we have found, easy to implement on a digital computer. 
It also lends itaeU to special-purpose computing hardware, such as 
digital or analog differential analysei^. In addition, by virtue of tiie 
form of the algorithm, the parameters that define the shapes are ex- 
tremely easy to compute. (In some cases they may require no com- 
putation at all. ) 

We intend to develop a method to construct complex arbitrary surfaces by 
piecing tegetiier surface "patohes. " Each such pateh will be defined by four 
boundary curves, in principle, althou^ it is harmless for one of the boundary 
curves to be degenerate, and to appear as a point instead of a curve segment. 
In the design of a surface, it is Intended that the designer begin with a single 
surface patoh, or a very small number of patehes, and then subdivide these 
regions with additional design curves defining boundaries of smaller patehes 
only when the internal surface needs modification. This is somewhat at variance 
with the customary procedure for mathematical curve fitting and surface fitting 
of existing curves and surfaces, in which a relatively large number of surface 
points already defined by some other procedure are used to obtain mathematical 
expressions for a surface that best fits them. Instead, the system to be de- 
scribed is intended to be used by the designer at the outeet, in the process of 
designing the surface, rather than later on as a means for making it mathe- 
matical. 

This is not to say that the surface-pateh technique cannot be used to for- 
mulate pateh-wise mathematical expressions for existing surfaces, but rather 
to indicate that the primary purpose of this surface technique is to facilitate 
the initial design process itself. 

When the design process is completed, the surface will be completely 
mathematically defined, since this definition occurs automatically and concur- 
rently witii design. 

Ordinarily a ship's hull or an airplane fuselage is described by certain 
important curves such as, in the case of the hull, a keel curve, a midships 
section, and a curve representing the sheer or deck line; these curves are 
sufficient to determine a surface, since they form the boundaries of a surface 
pateh. However, ordinarily this primary surface will not have certain desired 
characteristics, and it will have to be modified by introducing additional 
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information, such as for instance one or two other section curves. When these 
additional curves are introduced, the surface algorithm permits the computer, 
to "re-fair" ihe original hull form to contain these curves. 

Similarly, an airplane fuselage can be designed by drawing a profile cur\'e, 
a maximum half-breadth curve, and a mid-section of the fuselage. Again these 
curves suffice to define a primary surface, which in most cases will require 
modification by the addition of a few more curves to make more explicit tiie 
designer's wishes. As fliese additional curves are introduced, the original 
surface will be sub-divided Into patches, but the algorithm will automatically 
insure continuity of surface slope and curvature (if desired) and will incorporate 
these additional curves Into the surface automatically. This should make initial 
surface design virtually painless, and is intended to remove the tedious process 
of surface falrtng as It is now practiced in naval architecture. Airplane fuse- 
lages are usually somewhat simpler shapes than, say, yacht hulls, and for a 
nimiber of years second-degree curves have been successfully used for fuselage 
design; on the other hand, naval architects have steadily resisted the use of such 
methods In their work, since the complexity of yacht shapes makes it necessary 
to pay attention to the irksome details of the geometry involved, and second- 
degree curves prove to be cumbersome in such applications. 

The system that is described in tiiis report is intended to furnish the flexi- 
bility that second-degree curve techniques lack, and to remove almost entirely 
the need for the designer to be an analytical geometer. With this system imple- 
mented on a computer, there is reason to believe that the computer can take 
over all of the geometrical and mathematical burden of the design process, and 
leave the user free to be a sculptor assisted by an exquisitely skillful mechan- 
ical slave. 

Ultimately, when a graphical input-output hardware for a computer is 
available in the engineering design office, these methods will permit designers 
to delineate complex shapes with great ease, by simply drawing the salient 
curves that define and describe them. Already experiments along these lines 
are in progress in a few isolated laboratories both in universities and in 
industry. Very soon the two severe handicaps that have inhibited the wider 
use of such graphical devices will be removed. These inhibiting factors have 
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been hiii,h cost for the terminal hardware and small size of the working area. 
Rapid strides are being made on both these fronts, and within a few years it will 
be possible not only to draw on a virtually unlimited drawing surface , but to 
draw objects directly in three-dimensional space, and to view these constructetl 
objects as one would view an actual physical thing. 
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NOTATION 



We shall in what follows relate the x, y, and z coordinates of points on a 
surface to two independent variables u and w, so that we could write 

X = f (u, w) 
y = g(u, w) 
z = h^, w). 

If the functions f, g and h were specified, then for a pair of values of u and w, 
a point in space would be defined. If we held one of the independent variables 
fixed, say w, then by allowing u to vary, the point in space would trace out a 
curve. If subsequently we set w to a new fixed value and again allowed u to 
vary, we would trace out another curve, and so on. Clearly by stepping the 
values of w by small increments and allowing u to vary after each such step, 
we could produce a family of space curves that would lie on the surface and 
define it. All that is needed is some convenient and systematic way of arriving 
at the functions f, g, and h. 

It will turn out that the form of all of these three functions is the same; 
only certain internal numerical values are different. In vector notation we can 
write 

[x y z] = [f(u,w) g(u,w) h(u,w)] 

Sinoe^V^= fxyzj is a suitable conventional abbreviation for the vector 

on Ifae left, we Introduce a similar abbreviation for the right hand side: 

(uw) = [f (u,w) g (u,w) h (u,w)] 

Here, in the abbreviated symbol on the left, we shall omit the comma between 
the two letters. Later on, when no ambiguity can arise, we shall omit the 
parentheses as well, and write simply uw to stand for the vector. It is to be 
remembered that uw does not stand for the ordinary product of the two quan- 
tities, but is merely a bi-literal S3anbol standing for a vector whose components 
are functions of the two variables. 
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We plan to build up surfaces by adjoining surface "patches, " in an analogy 
of the piecewise fitting of complicated curves by curve segments suitably joined 
together. Accordingly, we shall at the beginning focus our attention on one such 
surface patch. To simplify arithmetic, we shall stipulate that the Independent 
variables, or parameters, u and w can take on only values between and 1. 
Then a surface patch can be considered to be a surface segment bounded by 
four space curves, (0 w), (1 w), (u 0) and (u 1). 




Here, typically, the symbol (Ow) stands for the vector describing the x, y, and 
z coordinates of points along the curve generated by allowing w to vary conti- 
nuously from to 1, while u Is held fixed and equal to 0. 

We shall introduce two scalar functions, F and F^ each a function of a 
single variable. These will be referred to as "blending functions" for reasons 
that will become clear. 

In order to compress the surface equation, and the proofs that we wish to 
demonstrate, we shall use a kind of indicial notation; we introduce the indices i 
and j, which can assume only the values and 1, and we invoke the customary 
summation convention for terms with repeated indices. This convention in our 
case simply means that when an index is repeated in a term, we write out all 
the possible terms that the actual indicial values generate, and then add them. 
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THE SURFACE EQUATION 

With these conventions and notational peculiarities in mind, we write 

(uw) = (iw)Fj<u) +(uj)F^(w) - (ij)F^(u)F.(w). 

(Typically, the first term on the ri^t expawis as follows: 

(iw)Fj{u) =(Ow)Fq(u) +(1w)F^(u). 

Thus the complete expansion would consist of eight terms, if carried out. ) We 
shall proceed to demonstrate that this surface equation represents a surface 
that contains the four boundary curves, and is thus defined by them. 

We must make a stipulation, a weak one, on the nature of the blending 
functions F and F : 

Fp(0) = 1 Fjj(l) = 



Fj(l) = 1 F^(0) = 



A further stipulation is that F and F be continuous and monotonic over the 
Interval. 

Now set u = a, where a can only be either or 1. Then, substituting in 
the surface equation, 

(aw) = (iw)Fj(a) + (aj)F.(w) - (ij)F^(a)F.(w). 

Consider F (a) vbidi occurs twice in the equation. By the stipulation, if i = a, 
F^(a)=l. 

Otherwise, if i ?! a, F.(a) = 0. 

Hence all terms in the expansion that contain i s^ a vanish; we can set i = a and 
what remains is 

(aw) = (aw)F (a) + (aj)F.(w) = (aj)F„(a)F,(w) 

= (aw) + (a})F.(w) - (aj)F.(w) 

= (aw). 
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This shows that for a = or 1, and hence (aw) = (Ow) or (Iw) , the surface equa- 
tion reduces to an identity. This implies that the surface contains its boundaries. 
An enttrely parallel argument would show that the equation also reduces to an 
identity for the other two boundaries (uO) and (ul). 

Provided a pair of functions F and F are chosen once and for all that 
satisfy tiie stipulations, Uie surface equation may be constructed immediately 
and uniquely for «iy set of boundary curves (uO) (ul) (Ow) and (Iw). It is to be 
observed Ifaat no restrictions have been placed on the form of the boundary 
curves; there is periiaps the restriction that they form a closed boundary, at 
least at the comers (ij) = (00), (01). (10), and (11) otherwise there will be mul- 
tiple values within the surface segment; similarly fliey should be continuous 
functions, but apart from these rather obvious restrictions, they can be of any 
shape whatever, including curves that can only be represented by tables of 
values . 

We can gain intuitive insight Into the nature of such a surface if we look 
at one of the terms, say (uj)F (w). 

We have the expansion 

(uj)F.(w) = (uO)Fq(w) + (ul)Fj^(w). 

This represents a weighted average of the quantities (uO) and (ul). When 
w = 0, F (0) = 1 and F (0) » 0, and the expression becomes simply (uO). As w 
increases, the wei^t of F^fw) decreases, while that of F^(w) increases, so that 
the surface partakes of Ae nature of both boundary curves. As w approaches 
the value 1, the influence of (uO) on the shape of the surface gradually disappears, 
while the influence of (ul) gradually becomes dominant. Finally, at w = 1, the 
curve (ul) represents the shape of the surface. We can say that the surface is 
generated by a gradual transition from (uO) to (ul), and that these two curve 
shapes are "blended" together by virtue of the blending functions F^^ and 7^. 
This discussion is some>»*at oversimplified, since we have omitted the term 
(iw)F. (u) and it too plays a part in determining the shape of the internal surface, 
as does of course the term involving the comer coordinates, (ii)Fj(u)Fj(w). 
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The entire surface equation is seen to be symmetric in u and w, and by 
virtue of this and a secondary symmetry in the functions F^ and F^ , we can 
abbreviate proofs about the behavior of the surface along all boundaries by 
exhibiting a typical proof for any one boundary, 

3 . 1 BOUNDARY S LOPE CONTINUITY 

It is our aim to design and delineate complicated surfaces by adjoining 
surface patches, in a pteeewise fashion. Consider two such patches A and B, 




with a common boundary. For patch A the boundary is (Iw); for patch B it is 
(Ow) , and the vectors of coordinates are equal , 

A (Iw) =B (Ow). 

Then the two patches will be continuous across their common boundary. They 
will however In general be discontinuous in slope across the boundary , and we 
wish to investigate this and make some amendments that will correct this dis- 
continuity of slope . 

We take the partial derivative with respect to u: Our symbolism for this 

partial derivative is (uw) = , '"^' , and when we substitute , say , u = , we can 
^ 'u u 

write (Ow) to mean the value of the partial derivative so obtained. Then 

(uw)^ = (iw) F| (u) + (uj)^Fj(w) - (ij) F| (u)F.(w) . 
Now substitute u = a = or 1 , as before . 

(aw)^ = (Iw) F| (a) + (aj)^F.(w) - (Ij) Fj (a)F.(w). 
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If we now place additional constraints on the blentUng functions , that their first 
derivatives 

F'(a) = (a = either or 1) 

we obtain the result 

(aw)^ = (aj)^F^(w), 
all other terms vanishing. 

This implies , for example , that when a = , 

(Ow)^ = (OO)^Fjj(w) + (OJO^F^(w), 

or , the derivative anjnuiiere along the boundary in the u direction (across the 
boundary) depends only upon the derivatives at the end-points of the boundary; 
it is entirely indepewlent of the shapes of the four boundary curves , including 
the boundary (Ow) itself. '^ 

Thus for the two patches A and B, if 



and 



A(10)^ = B(00)^ 



A(ll)^ = 3(01)^ 



1. e. , if the boundary curves are continuous In slope in the u direction at the ends 
of the contiguous boundary between patches , we are guaranteed to have 

A(lw) = B(Ow) everywhere along the boundary regardless of the 
shapes of the boundary curves of A and B. This is a remarkably powerful and 
useful property , achieved at the sli^t expense of extending the stipulations on 

the F.. 

1 

Similarly, the secc»id derivative with respect to u is 

(uw)^^ = (iw) Fj- (u) + (uj)^^F^(w) - (ij) F'' (u)Fj(w) 

and if we further stipulate that F|' (a) = we obtain 

(aw) =(aj) F.(w). 
* 'uu uu i 
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This establishes second derivative (or curvature) continuity as an auto- 
matic and inherent property of adjacent patches, provided their boundary curves 
have this kind of continuity at the end-points of the boundary. It is easy to see 
that we may escalate in this way to any level of derivative contlmiity we wish 
along contiguous boundaries . 

3.2 SLOPE CORRECTION SURFACE 

The surface equation already described is very general , in the sense that 
it can contain virtually any boundary curve we wish, and it has certain benign 
properties of derivative matching along boundaries ; nevertheless it is not a 
universal formula for all surfaces , and there are many that do not belong to its 
family. We have already seen that surfaces generated by the surface equation 
have a definite intrinsic slope along boimdaries , «^ose variation is rigidly 
prescribed by a simple formula. Obviously surfaces exist whose boundary 
slopes do not match this Intrinsic slope , except at the end-points of boundaries . 
Nevertheless , we wish to be able to patch together such other surfaces with our 
special surfaces , so as to have slope continuity (or continuity of any level of 
derivative). 

To do so, we Introduce a new surface equation, describing a slope-correc- 
tion surface , which when added to the first surface equation has the property of 
leaving the boundaries unchanged , but causing the derivatives across boundaries 
to vary in any arbitrary way we wish, as we move along the boundary. 

The equation resembles the first form very stron^y . It is 

(uw) =(iw)^Gj(u) +(uj)^G.(w) - (ij)^^G.(u)G.(w). 

Here, typically, (iw) is a function of w only, and describes the arbitrary 
variation of the derivative with respect to u as w varies , along the curve (iw) , 
and similarly for the other boundaries . The vector (Ij) represents the cross 



derivatives of the four corners. 


Typically , 


(00) = -iV^ 

'uw ouow 


u =0 
w = 
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The functions G and G are again blending functions or weighting flinc- 
ttons, but they have properties different from the ftinctions F^ and F^^. We 
stipulate 

G(,(0) = G^(0) = Gq(1) = G^(l) = 

GJ|(0) = 1 G^(l) = 1 

GJj(l) = G^(0) = 

or in the indicia! notation used earlier, 
G (a) = , a and i = or 1 . 

GJ(a) - 0, a^l. 

G|(a) = l, a = i. 

We need to ensure that the vectors describing the boundaries vanish 
identically, and that the vectors describing the slope variation along boundaries 
are indeed given by the equation. The proof proceeds along precisely the same 
lines we used brfore. First, substitute u = a. The equation becomes 

(aw) = (iw)^Gj(a) + (aj)^G^(w) - (lj)^^Gj(a)Gj(w). 

Consider (aj) . We wish to have the correction surface leave the original 
boundary vectors unchanged , and hence the boundary vectors of the correction 
surface must vanish; ie, 

(iw) = 
(uj) = 
Then the derivatives of these boundaries must also vanish; in particular, 
(iw) = and then (ij) = , when w = j . Hence (aj)^ = . 

Thus (aw) = indicates the desired behavior of the correction surface along a 
boundary. Similarly (ua) = 0. 
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To examine the slope variation along a boundary , differentiate the 
equation with respect to u; 

(uw)^ = (iw)^GJ(u) + (uj)^^Gj(w) 

-(ij) g:(u)G,(w). 
'uw i * ' j 

Now substitute u = a 

(aw)^ = (iw)^Gj (a) + (aJ)„^G.(w) - (IJ)^^G! (a)G.(w) 

As before , G I (a) = 1 if and only if a = i , 
so we get 

(aw) =(aw) G' (a) + (aj) G,(w) - (aj) G' (a)G.(w) 
^ 'u 'u a ^ ' ^ ■"uw j ' uw a ' y 

= (aw) + (aj) G.(w) - (aj) G.(w) 
'u uw j uw ] ' 

= (aw) . 
' 'u 

This demonstrates that the surface has the slope variation along the 
boundaiy as required. To make use of this slope correction surface, we must 
first determine what the intrinsic slope of the surface to be corrected is , and 
then we must subtract this elope from the desired boundary slope , to yield the 
correction slopes that enter into the equation. Thus if (Ow) is the desired 
slope, and 

I (Ow) is the intrinsic slope , then 

C (Ow) will be the correction slope , 

C (Ow)^ = (Ow)^ - I (Ow)^. 

The correction slopes C (iw) and C (uj) are the four functions that enter 
into the slope correction surface . . The desired surface is obtained by adding 
the correction surface to the first surface : 

(uw) = I (uw) + C (uw) 

where we use the symbol I (uw) to represent the surface whose boundary slope 
is being modified. 
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3,3 HIGHER-ORDER CORRECTION SURFACES 

Analogous forms may be obtained for correction of hi^er derivatives along 
boundaries. For second derivative correction, the surface equation is 

(uw) = (iw)^Hj(u) . (uj)^Hj(w) - (iJ)„^^Hj(u)H3(w). 

In this equation, the blending ftmctlons H have the stipulations that, for a = 
or 1 as before, 

Hj(a) = 

H.'(a)=0 

Hj"(a)=0,i?^a 

H^"(i) =1, i=a. 

With these constraints on the H , it is easy to arrange matters so that this 
second-order correction surface is zero everywhere on the boundary , has zero 
slopes across boundaries, and has second derivatives across boundaries speci- 
fied by (iw) and (uj) v^atever these functions may be . The addition of this 
surface vector to a given surface vector will then provide a means for boundary 
second-derivative correction without disturbing either the boundary shapes or 
boundary slopes. 

Althou^ we have already carried out a similar proof for slope correction, 
it mi^t be well to exhibit once again the course of the argument. 

First, to show that the boundary vectors are zero, substitute u = a: 
(a = or 1.) 

(aw) = (iw)^^Hj(a) . (aj)^Hj(w) - (iJ)„„^Hj(a)Hj(w) 
= (aj)^Hj(w). 

The term (ai) refers to the second derivative in the w sense at each of the 
* "ww 

four corners, such as, typically, (00). As in the case of slope correction, we 
must have 

(iw) = along boundaries . 
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Then (Iw) = 0, (ij) =0, and in particular (aj) = 0, so that the 

•nrf WW " WF 

equation satisfies the boundary condition. 

For boundary slope vectors, differentiate with respect to u: 

I t 

(uw) = (iw) H (u) + (uj) H,(w) - (ij) H. (u)H.(w) 

'u ^ 'uu i ^ ' ^ -"wwu j* ' ^ "uuww 1 * ' y ' 

Set u = a: 

(aw) =(aj) H.(w). 
'u * •"wwu y ' 

We wish to have the slope vectors vanish along boundaries , so typically 

(Iw) * for all w. 
'u 

But then (Iw) = and (iw) = by taWiw derivatives . The order of 
uw * 'uww 

differentiation is immaterial , so 

(iw) = (iw) , and finally we can conclude that 
'uww 'wwu 

(aj) = 0; again the ri(^t and left hand sides of the equation are 

in agreement . 

Finally, we differentiate again with respect to u: 

(uw) =(iw) H (u)+(uj) H,(w) - (ij) H, (u)H,(w). 

^ 'uu ^ 'uu i ^ ' * ■"wwuu j' ' ' -"wwuu i * ' j^ ' 

Set u = a; only terms in «4Uch a = i remain: 

(aw) =(aw) H"(a)+(aj) H.(w) - (aj) H (a)H.(w) 

* 'uu 'uu a 'wwuu J ' •^'wwuu a ' j ' 

= (aw) 
'uu 

Again we have demonstrated an identity . The escalation to any level of 

boundary derivative correction vector is obvious . 

3.4 MATRIX FORM 

The surface equation 

(UW) = (iw)F.(u) +(uj)F.(w) - (ij)Fj(u)F.(w) 
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may be expanded directly into matrices, to yield: 



(uw) 



[uO ul] 






[f^u F^u] 



Ow 
Iw 



- [^o'^ ^A 



00 
10 



01 

11 



F^w 
Fj^w 



In this we have treated the Indicial form term by term in a straight- 
forward way. We shall in what follows omit parentheses, since no misunder- 
standing can arise. Thus typically F u is written in place of F^{n) as a matter 
of convenience and economy. Similarly, typically 00 is written instead of (00); 
the reader should be reminded that this is merely a compact way of exhibiting 
the X , y , z coordinates at point (00) . 

It means: 

00= rx(OO), y(00), z(00)] when written out completely . 

The three vector (matrix) products are equivalent to the following three products: 



[l F^u F^u] 



[l FqU F^u] 



+ [l FpU F^u] 



uO 







Ow 

Iw 



ul 




1 







F^w 







F^w 








-00 
-10 



1 











1 







F^w 







F^w 




-01 

-11 



FqW 
Fj^w 



and in this form we can perform the addition, obtaining 



(uw) = [l FgU F^u] 






uO 


Ow 


-00 


Iw 


-10 



ul 




1 


-01 




F„w 


-11 




M 
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It is subtly more convenient to rewrite this in the equivalent form 



(uw) = - [-1 FqU F^u] 






uO 


Ow 


00 


Iw 


10 



ul 




-1 


01 




F^w 


11 




. V. 



so as to avoid the awkward minus signs in the 3x3 matrix . 

Two facts should be noted. The leading row vector in front of the matrix 
and the trailing column vector following the matrix are transposes of one 
another , but with different arguments ; the matrix represents the boundary con- 



ditions of a patch , The partition 



is redundant, since its elements 



00 01 
10 11 
must agree with uj and iw for u and w equal to or 1 . 

We have already suggested that we can maintain slope continuity across 
boundaries by suitable stipulations on F , and we have also already suggested 
that ^en desired we can adjust slopes across boundaries by a second additive 
vector with suitable stipulations on its G. . We shall now investigate the com- 
bined form of the surface equation. To do so we shall prefix a symbol to the 
vector uw to indicate whether we are talking about the first surface equation, 
or the coirrection surface equation, and we shall omit the prefix symbol when 
we are talking about the combined form. Thus 

uw = suw + cuw , with 
suw = the primary surface 
cuw = the correction surface 

uw = the combination. 

Accordingly, using this notational convention, we will take derivatives, with 
respect to u , of the surface equation suw in order to determine its slope vector 
in the u direction. 



20 



SECTION m 



suw = ■ 
u 



[O Fg'u f/u] 



[-1 F^u F^u] 



We substitute u = , and obtain 

sOw = [sOO 801„| 
u L u OJ 



suO sul 

sOw sOO sOl 
slw slO sll 



-1 
F^jW 
F w 











F„w 



suO sul 

u u 





-1 




FqW 




LVJ 



Fj^w 



Now consider, for example, sOO and the desired 00^. The symbol sOO^ 

refers to the slope vector at a corner; we have already seen that at corners the 

correction surface cOO = 0, and so 00 = sOO . This is bourne out Intuitively 
u u u 

by the reflection that at (ij) comers , the two crossing boundary curves com- 
pletely define the slopes there; since this is so, no correction of slope need or 
can be applied . 



Hence we should write 

sOw = foo 
u i u 



-J 



F w 



By analogy and symmetry we can write the remaining three statements ; 



slw 



BUO 



sul 



[lO 

I- u 

[v 



uJ 
F^u] 

F^u] 



F^w 

I F w 



00 

V 

10 

L w J 



01 
y 

11 



/: 
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In order to obtain a desired slope vector along any of the boundaries , we add 
the correction surface, whose equation is 



-[- 



1 GqU G^u 



cOw 
u 

clw 



cOO 

u 

clO 



uw 



cul 
cOl 



ell 



L 
w 




-1 


L 
UW 




G^w 


L 

UW_ 




G^wJ 



As we have already remarked , the desired surface uw is the sum of the vectors 
suw and ouw. Hence the correction slope vector, such as, typically, cuO , is 



cuO 



= uO - rF„u F,ul 
w L 1 J 



00 

V 

10 



This is an entry in the correction surface matrix. 

Now we introduce a new fact: the comer cross derivatives of the primary 
surface equation are aU zero. To show this, differentiate the indlcial expres- 
sion first with respect to u, then with respect to w, and finally set u = a, w = b, 
where a and b are as usual either or 1 . We have 

uw^ = (iw)F^ (u) + (uj)^F^(w) - (ii)F^ (u)Fj(w), and 



(iw)^Fj' (u)- + (uj)^Fj (w) - (ij)Fj' (u)F.' (w). 



Evidently this expression vanishes for (uw) = (a b) . This shows that the corner 
"twists" , or cross derivatives , of the original surface all vanish; it is a 
peculiarity of the first fundamental surface equation. 

Hence we can assert that 

cij = ij ; This says that the desired twists at corners are 
uw uw 

identical with the correction surface twists , since the fundaments surface has 



no twist . We shall use this result to replace the partition 



cOO 



clO 



cOl 



ell 



with 



00 



01 



uw 



10 



11 
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We can rewrite the expression for cuo as follows: 

w 



cuo 



[1 -F^^u -F^uJ 



and 



cul = 



1 -F^u -F^uJ 



and of course, perhaps trivially, 



= 



1 -F^u -F^ul 



uO 

V 

00 

V 

10 



ul 

V 

01 

11 









Each of these matric products represents an element of the top row of the 
correction surface matrix. 

Since the row matrix ^X -F u -F u] is common to these three products, 
it can be factored out and introduced into the matrix |^-1 G u G uj to yield 

[- V -*^o" -^1^ ^'0^ ^1^ 
which is the same as the vector 

[-1 F^u F^u G^u G^u] 

We replace the elements of the top row of the correction surface matrix 



by the three matrices 








'uO 
w 







00 
w 







10 
L wj 



and 



ul 
\ 

01 
\ 

11 
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This causes it to become a 5 x 3 matrix, and we now have the intermediate 
result, 



-1 F^u F^u G^u G^u 






uO 
w 


ul 

w 





00 
w 


01 

w 





10 


11 

w 


cOw 
u 


00 
uw 


01 
uw 


clw 
u 


10 
uw 


11 

uw 



-1 



G^w 

G w I 

. ^ J 



By similar procedures, we can write for the elements of the first column 
of the correction surface matrix, 



cOw = Tow 00 01 1 
u L u u u J 



clw 



1 



•F^w 
■F w 



= fiw 10 11 1 r 1 

u L u u u J 



-F^w 



-F w 



and again trivially , perhaps , 



^[ 



0=000 



-V 



^F w 



When we factor out the common column matrix as before, and replace each 



entry of the column matrix 









cOw 
u 

clw 
u 



by the above expressions, 
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we obtain the complete matrlc expression for the correction surface: 



cuw = - [-1 F u F^u G^u Gj^uJ 












uO 
w 


ul 
w 




-1 











00 
w 


01 
w 




V 











10 
w 


11 

w 




F^w 


Ow 
u 


00 
u 


'\ 


00 
uw 


01 
uw 


Iw 
u 


10 
u 


11 
u 


10 
uw 


11 

uw_ 




.V 



If now we border the original surface equation matrix, it can be written. 



suw = - [-1 F u F u G u G u ] 






uO 


ul 










-1 


Ow 


00 


01 










V 


Iw 


10 


11 










F w 
G„w 


































G^w_ 



In this borderi:^ process, the value of the matric product is unchanged. 

Since the pre- and post-multiplicative matrices in this equation are the 
same as those of the correction surface equation, we can add the two 5x5 
matrices and pre- and post-multiply by the two vectors. We shall perform, in 
fact, 

uw » suw + cuw, and obtain 

uw = - [-1 FqU FjU GjjU Gj^u j 






uO 


ul 


uO 
w 


ul 
w 


Ow 


00 


01 


00 
w 


01 
w 


Iw 


10 


11 


10 
w 


11 

w 


Ow 
u 


00 
u 


01 
u 


00 
uw 


01 
uw 


Iw 

u 


10 
u 


11 
u 


10 
uw 


11 

uw. 





' -1 




FqW 




F w 




G^w 




.V 





This is a general expression for a slope-matching, slope contimious surface patch 
with entirely arbitrary boundaries and entirely arbitrary slopes across these 
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boundaries. There are no stipulations whatever on the nature of the boundary 
slope function. The stipulation on the F and G functions have already been dis- 
cussed. 

Now that we have constructively arrived at a general expression for sur- 
faces that have a prescribed boundary vector and a prescribed boundary slope 
vector, it mi^t be interesting to apply a proof to a conjectured higher order 
surface equation in which not only boundaries , boundary slopes , twt also 
. boundary second derivatives are vector quantities under control . 

We postulate, therefore, that by analogy the surface equation is 

uw = - f-l F u F u G u G u H u H u] 

T r 
-1 

F w 
X Ow 00 01 00 01 00 01 G w 

G w 
HqW 
H^w 

It represents a surface patch whose vectors of coordinates , slope , and curva- 
ture as well , are everywhere arbitrary along its boundaries . The first colunm 
ami first row of the 7 x 7 matrix represent these boundary conditions; the re- 
mainder of the matrix is redundant, since the quantities this partition contains 
must all come from the colunm and row by differentiation. 

We can test this equation by seeing whether it contains a boundary curve . 
To this end, set u = 0, so that we check whether it contains the boundary (Ow). 
We obtain, invoking the stipulations on the F, G, and H functions. 






uO 


ul 


uO 
w 


ul 
w 


uO 

WW 


ul 

WW 


Ow 


00 


01 


00 
w 


01 
w 


00 
uww 


01 

WW 


Iw 

Ow 

u 


10 


11 


10 
w 


11 

w 


10 

WW 


11 

WW 


00 
u 


01 
u 


00 
uw 


01 
uw 


00 
uww 


01 
uww 


Iw 
u 

Ow 
uu 


10 
u 


11 

u 


10 
uw 


11 

uw 


10 
uww 


11 

uww 


00 
uu 


01 
uu 


00 
uuw 


01 
uuw 


00 
uuww 


01 
uuww 


Iw 
uu 


10 
uu 


11 

uu 


10 
uuw 


11 

uuw 


10 
uuww 


11 

uuww 
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Ow = 



[- 



110 



00 01 00 01 

w \ 

Ow 00 01 00 01 



00 



00 



01 



01 



F^w 



F w 



G^w 



G w 



H^w 



H w 



In the boundary matrix we have omitted irrelevant terms, because of the 
zero's in the pre-multipljlng vector. We obtain, by performing the multipli- 
cation, 

Ow = - [ow o] 



-1 



F w 




F w 



G^w 



G^w 



H^w 



H w 



-'= Ow, which is the hoped-for identity. 
We can next try to see whether the equation also conforms to the boundary 
second derivative conditions . It will be more convenient in what follows to 
introduce some abbreviated notation. 

Set 



fu 



and a similar expression for [fw 



.1 F^u F^u GgU G^u H^u H^u 
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Set the 7 X 7 boundary condition matrix equal to I Bj . 
With these abbreviations , the surface equation is 

uw = - [fu] [b] [fw] '^. 
We differentiate with respect to u: 

'-u=- (H M ^ H W) W 



and again: 



uw 



--(H N -2 [fu] [bJ . [fu] [Bj)[fw] 



We wish to investigate the rig^t hand side of this equation for u = , that is , 

for uw = Ow . The blending function vector and its derivatives become 
uu uu " 



If'O 

If'O 



= [0000016 

= [0001006 

[fo] =[-1100006 

As for the first and second partial derivatives of the [bJ matrix , all 
elements of [b ] and [b 1 vanish except for those in the top row . 



Then 



[f'ol fBl = fow 00 01 

L J L J L uu uu 1 

2 [f'o] [b J « [ the null vector] 

[fo] [b ] =[ -00 -01 

L J I. uu-1 I- uu 



00 01 00 01 

uu uu uu uuw uuw uuww uuww 



-00 -01 -00 -01 

uu uu uuw uuw uuww uuwwJ 



The sum of these vectors is evidently 



Ow 

'• uu 



Finally , 



Ow 



-I' 



Ow 



0] 
000000] [fw] 



Ow , as expected . 
uu 



28 



SECTION m 



We have shown that the extended surface equation satisfies the second 
derivative boundary conditions . In a sinallar way it can be shown to satisfy the 
first derivative boundary conditions, but this was skipped In favor of the proof 
for the hl^er derivative , since the procedure exhibits a few interesting points , 

By analogy we could construct matrix products to represent surfaces 
which satisfy even hi^er derivative conditions across boundaries . 

3.5 BOUNDARY CURVES 

It is often convenient to use particular boundary curve functions defined 
by the curve end-points and end-point tangent vectors . We can use the blending 
functions themselves to define such curves . For example , the uO boundary 
curve can be described by the equation 



uO = [FqU F^u GjjU G^u] 



00 

10 

00 
u 

10 



\>diere the colunon vector contains the end-point Information. We observe that 
the row vector becomes [l o] vtlien u = 0; it becomes [o 1 oj 
when u = 1. Again, if we take derivatives of this row vector with respect to u 
we obtain Ff 'u F "u GL'u G 'u] and this becomes [o 1 oJ foru = 0, 
and it becomes lo ij foru = l. 

With this behavior of the row vector , it is easy to see that the equation 
does indeed represent a curve satisfying the end-point conditions . 



The matrix form of the surface equation has been shown to be 



uw = - [-1 F u F u G u G uj 






uO 


ul 


uO 
w 


ul 
w 


Ow 


00 


01 


00 
w 


01 
w 


Ow 
u 


00 
u 


01 
u 


00 
uw 


01 
uw 



Iw 10 11 10 11„„ 

u u u uw uw 



-1 



F^w 
G^w 
G w 



<^ 
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Now when, in computing uw, we perform the matrix multiplication from 
the left, we have for the second column 



[-' ^" ^" S^ ^1^ :^^° 

00 
10 
00 



10 



-uO 



'■'O^^ FjU G^u G^ul 



00 
10 
00 



-uli - uO ^ 0. 



10 



Now if similarly ul, uO and ul are functions of the same kind, their 
w w 

corresponding column products vanish just as in the case of uO. .4ecordini;iy , 
the resulting product of the three matrices has the form 

r 







-1 



Y w 





F w 



'^^o^^- 



.^^1^^: 



where P is the product of the row vector and the first column of the malri.\, or 



1 F^^u F^u G^^u GjU, 



1 



Ow 
Iw 
Ow 



Iw 



F^u F^u G^u G^uJ 



I Ow 

llw 

Ow 
u 

Iw 



/: 
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If again the elements of the column vector, Ow, Iw, Ow^ and Iw^ are 
likewise ftinctions described as outlined, we can write typically 



Ow 



[ 00 01 00 01 ] 
' w w' 



V" 


F w 


G„w 


LvJ 



For the complete column vector we have 

Ow = 00 01 00.. 

Iw 



Ow 
u 



iw 
. u_ 



00 


01 


10 


11 


00 


01 



10 



11 



10 



00 



10 



01 



11 



01 



11 



w 




V 


w 




F w 


uw 




V 


uw. 




G^w 



When we substitute this result for the column vector, we obtain the surface 
equation 



uw = [ F.u F u G u G u ] 



00 


01 


00 
w 


01 
w 


10 


11 


10 
w 


11 

w 


00 
u 


01 
u 


00 

uw 


01 

u\ 


10 

L u 


11 
u 


10 
uw 


11 

u\ 



F w 



G^w 
G w 



This is a particularly convenient form for computation. The 4x4 
matrix contains nothing but information about the corner coordinates, corner 
slopes, and corner twists; all entries are constants, and the partitions of the 
matrix systematically group these quantities. The leading row vector and 
the trailing column vector are transposes of one another, (but with different 
arguments, of course. ) 

We shall refer to the 4 x 4 matrix as the "boundary condition" matrix, 
and shall assign to it the symbol B, so that the matric equation for the surface 
could be written 



/; 
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uw = [ FqU F u G u G u] B 



r~ "1 


F„w 


F^w 


GqW 


k-J 



It must be remembered that each of the entries in B is a three-vector, 
whose components are x, y, and z coordinates and slopes and twists. This 
means that B is really a tensor. 

3.6 BLENDING FUNCTIONS 

We can relate the blending function vector to a so-called basis vector 
in the following way. Let { u u„ u„ u ] be a vector whose elements are a set 

X ^ o 4 

of linearly independent functions of the variable u. Then we can postulate 
the existence of a matrix M such that 



fV V V S"J = f"i "2 "a "4] 



M. 



To evaluate the M matrix, we substitute u = 0, u = 1 on both left and right 
hand sides of the equation. Then we take derivatives of both sides, and 
again substitute u = and u = 1. There results 



V ^1° V °i°" 


= 


'\ 


°2 S 


< 




F„l F^l G^l G^l 




h 


h h 


h 


M 


FJjO F^O GJjO G'O 




'[ 


°2 °3 


°; 




f;,1 F'l G'l G'l_ 




}l 


'2 ^3 


\ 




matrix on the left is the identity matrix 


10 





by virtue 




1 









1 
















1 
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of the stipulations on the blending functions and their derivatives. From this, 



we 


conclude that 






^0 


n 


n 













"l 


2 


3 


4 










1 


1 


1.. 


1 , 










1 


'f. 


■A 


4 








-1 














M 




0' 


n' 


0!. 


"', 










1 


2 


3 


4 










T 


V 


i:. 


1' 










. 1 


2 


3 


4 



and we need only find the inverse of tiWs matrix (if possible) to obtain M. 



(In the matrix, the notation 1' means A(^^ 



du 



typically. ) 



u=l 



In the next section we shall for the first time be specific about the 
basis vector [ u u u u ] , but it is interesting and important to realize 
that so far in the discussion nothing has been said to diminish the generaUty 
of the mathematical structure. It Is hoped that ttie reader will not lose 
sight of the fact that the surface equations in their several forms can be 
implemented in many ways. We propose to develop one such implementation 
in detail, but it is only one of many. 

3. 7 CUBIC BASIS VECTOR 

Let the basis vector be 

f"l "2 "3 "4^ " f" " " ^'' 
The vector on the right contains four specifically chosen linearly independent 
ftmctions of u. the powers of u, and when multiplied by a coefficient vector 
yields cubic polynomials: 



,32 

[ u u u 1 ] 



Au^ + Bu^ + Cu + D, 
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By the reasoning of the last section, we have for this vector basis, 



M-l = 



1 

1111 

10 
3 2 10 

from which we can obtain the desired inverse 

2-211 

M = -3 3-2-1 

10 

10 



Now we can write 

3 2 
[ FjjU F u GqU G u ] = [ u u u 1 ] M. 

We shall abbreviate the notation for the basis vector in what follows. We 

shall write 

3 2 
[u u u 1] = U 

3 2 
and [ w w w 1 ] .= W. 

The matrice surface equation 

uw = [ F u F u G u G u ] B 



'V" 


F w 


V 


G w 

. 1 



now becomes, simply and compactly, 
uw = U M B M^ w\ 



(Superscript t means transpose. ) 



If U and W are cubic basis vectors, then the surface patch is the so-called 
bi-cubic surface. Such surfaces are very easy to compute, particularly 
since the basis vector is so easy to evaluate. In passing it is important to 
remark that the above compact surface equation is not limited to cubics; 



yr 
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U and W are not restricted to cubic basis vectors, and M is simply the matrix 
that generates the appropriate set of blending functions. Among other pos- 
sibilities, U and W might be higher order polynomial basis vectors; or they 
might be any set of linearly independent functions. Provided the associated 
m"'^ matrix has an inverse, these basis vectors are acceptable. 

We can write, for w held fixed, an expression for a u - varying cvrve 
on the surface: 

uw = U (M B M* W*) = U A 

where A is a column vector of constant coefficients. We can write a simi- 

t 
lar exiH-essIon for u held fixed and w varyii«. The matrix product MBM 

is the same in either case. This suggests that for any surface patch this 

product should be evaluated first; thereafter, we can either obtain u-varying, 

w-constant curves or w-varying, u-constant curves in an obvious way. 

We shall investigate another basis vector that is composed of another 
set of Unearly Independent functions (not powers of u) in a later article. 

3. 8 DIFFERENCE EQUATIONS 

If the basis vectors are polynomial bases, we can invoke the techniques 
of finite differences to calculate points on the surface patch. 



Consider the matrix 



Then 



10 
110 
111 
111 

a 
b 
c 
d 



a 
a -)- b 
a + b -f c 
a + b + c + d 



If a, b, c, d are respectively third, second, first and zero-order differences 
of the cubic n', then the column matrix on the right of the equation represents 
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the corresponding differences for the cubic (n + 1) . The differences for 



(n + k) are given by 



udj 



where L means k successive 



multiplications by the L matrix. 

When n = 0, we can easily find that 



for cubic s. 



a 




^ a' 


b 




-6 


c : 


1 


. d. 




L0_ 



Using this, 



k = [0001] 



In this expression, the vector [0001] serves to select the bottom element 
of the resulting column vector after k multiplications by L. 



By extension, we can write the more complete statement 

6 



[ k""^ k" k 1 ] = [ 1 ] L 



-6200 
1-110 
1 

We shall call the [ 1 ] vector 1; and we shall call the 4 x 4 matrix N, 
so that 

r k^ k^ k 1 ] = 1 L*^ N. 
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Now let the usual parametric variable u be represented by u = kS , where 
k = 0, 1, 2 ... . and where 8 la an Increment size. Then 



3 2 3 2 

[u u ul]=[kk kl] 



S3 



82 



8 



Call this last square matrix A . Then 

u = il'^na. 

This expression states that we may step along the u parameter, in 8 incre- 
ments, by successive multiplications by the L matrix, and thus evaluate 
the U vector at these steps. In order for u to go from to 1, k must start 
at and go to -s- , since k 8 = 1. 



We can also write, for the W parameter 
W 
The sui^ace equation in difference form can now be written out in 



IL^NA 



full: 



IL^NA M B M* A* N* L* 1*. 



Call the partial product N A M B M A N = S... a square 4 x 4 matrix. 



Then 



and 



L S_ = S . a new square matrix. 



L'' S-Q = S , after j multiplications. 



We remember that for any column of S - the multiplication by L is a process 
of cumulative addition, as shown by 



a 

a + b 
a + b + c 
a + b + c + d 
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We can write in general that 

L S., L* = S. , , , 
]k 3 + 1, k + 1 

where the new square matrix is obtained from the old by cumulative addition 
of column elements, followed by cumulative addition of row elements. These 
operatlona are ftirthermore commutative, which means that we obtain the same 
result If we first add row elements and then afterward add column elements: 

We have, finally, that at u = 0, w = 
the siurface equation Is 

00 = 1 S-jj 1 ami in general 

uw = 1 S . 1 *^ where u = iS , w = kS , and S., has been formed 
jk JK 

from S by j column additions and k row additions. 

This obviously furnishes an extremely simple way to generate discrete points 

on a surface patch. The pre-multiplier 1 has the effect of selecting the bottom 

row of S.. and similarly the post-multiplier 1 has the effect of selecting the 

jk 
last column of S., . The bottom right hand corner element of S , is the value 

of the coordinate for a point on the surface, atu=jo,w = kS. 

Consider 1 S. . This represents the row vector obtained after j cumula- 
tive addition operations have been performed on the columns of 8^^^. The 
right hand element of this vector is the value of the coordinate at u = jS, 
w = 0. We can hold u fixed and step out successive values of the coordinate 
for w - varying, simply by cumulative addition on this row vector alone. In 
this case, the resulting right hand element is the marching coordinate value. 

T 
An analogous remark can be made for the product S^j^ 1 . This is a 

column vector, and successive cumulative additions of Its elements marches 

out values of the surface coordinate for w = kS fixed, and u varying. 

Althoi^h the arithmetic of the foregoing difference method is very 
attractive. It possesses certain drawbacks that must be made explicit. The 
coordinate values are precise if and only if no truncation error whatever is 
allowed in the arithmetic. Error is cumulative, and the least departure from 
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a precise number at the start will rapidly propagate. In the A matrix there 
are numbers of S magnitude; for 100 calculated points along a curve, this 
calls for 6 decimal digits, all exact, or 18 binary bits. For display purposes, 
100 points along a curve are adequate, but for engineering purposes this 
information is too sparse. To calculate intermediate points, the A matrices 
need to be changed and the S matrix recalculated. Furthermore, it is clear 
that we soon reach an upper limit on the number of available bits in the 

4 

computer word, becasue of the r^id growth of S . 

We can of course use the difference technique to calculate coefficients 
for a cubic and then calculate points using them. For consider 

uw = U M B M*^ A*^ n"^ l'^'' 1^. 

T T T Tic 
Jn this, the partial product M B M A N L represents a square 

matrix whose last column consists of the coefficients for the cubic at fixed 

w £: k 8 , with u varying: If S. is this square matrix, 

A 



uw = U S. 1 ^ = U 



B 
C 
D 

The cubic can be generated by digital methods for any u, or it could be 
generated quite easily by anak)g differential analyses hardware. The inte- 
grators of tbe analog differential analyses are loaded £rom the values of the 
cohinm coefficient vector, and it is then undamped and allowed to generate 
the curve. Meanwhile the digital machine can perform a cumulative addition 
on the rows of S. to yield S. . The last column of this new matrix then 
contains the new coefficients for the curve ofw= (k-t-l)S, ready to be 
loaded into the integrators. 

Recent developments in hardware, particularly a hybrid digital- analog 
multiplying device, may make it possible to generate surfaces for display 
directly £rom the iodiclal forms of the equation, or from the matri:£ form 

uw a U M B M^ W^ 
without the necessity for recourse to this last difference technique. K is 
however interesting to know that the method exists. 
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SECTION IV 

HYPERSURFACES - HIGHER DIMENSIONS 

We can readily extend the surface equation to describe hyper- surfaces 
immersed in hyper-space. For this purpose we shall introduce a slight variant 
on our notation. We shall write typically 

. = F.u. 1 = or 1. 

1 X 

This will be a standard replacement for the blending function notation. The 
stipulations on the F. are as before, so that if u = a, a = or 1, we can write 

. = when a / i 

. = 1 when a = i. 
1 

For slope continuity across boundaries, typically 

a' 
i 
function with the argument = a 



where this symbol means the first derivative of the blending 



For higher order continuity across boundaries, the additional stipulations 
on the blending fucntions are the same as for ordinary surfaces with two degrees 
of freedom and have already been discussed. 

The general surface equation for hyper-space is, in indicial form, 

(uvw . . . ) = (ujk ...)... 

^(iv....)^:... 
u v 



(ijw . . . ) . . 
1 } 



uvw 
-(N-l)(ijk...). .^... 



/ 
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In this equation, N is the number of independent parameters in (u vw. . . ); it is 
the number of degrees of freedom of a point on the hyper-surface. The indices 
i, j, k etc. can take on only the values or 1. 

Let us proceed to prove that this surface contains a boundary, say for 
example the boundary (u . . . ). We hope that the following equality holds: 

V 
(u . . .) ={ui k . . .). 1^ . . . 

, u 
+ (i k ...).,.. . 
1 k 

u 
- (i j . . . ) 



1 1 



u 
-(N-l)(i]k. . .)ijk- ■ • 

The last term in this expression is non-vanishing if and only if all indices other 
than i are zero, i. e. , j =■ 0, k = 0, etc. We can accordingly rewrite this term 
as 

(N-1) (i 0. . .)". 

Next consider the second term on the right; 
(iOk. ..)-... 

It is non-vanishing if and only if k = 0, etc. 
We can accordingly rewrite it as 

(iO . . .)|' . 

A similar consideration applies to 

\ (i i ...).. . which also becomes 

\ i J 

(i . . . ) I" . 
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There are evidently N-1 such terms, all identical, and they are removed by the 
last term. All that is left is the term 

but since J = and k = 0, we finally have (u . . . ) on the right. This establishes 
the identify, and the surface equation has thus been shown to contain this boundary 
curve. It is trivial to show that the surface contains all boundaries, and is defined 
by them. 

We can also show that the hyper-surface contains boundary surfaces 
of lower order. We shall content ourselves with the case for N = 3, and show 
that it contains surfaces for N = 2 which are identical with our ordinary surfaces. 
We have 

(u V w) = (u j k) J" ^ 

Mivk,;'- 

Set V = 0. Then substituting, and retaining only non-vanishing terms, 

V 
(which means that j must be replaced by whenever it occurs, and . = = 1) 

(u w) = (u 1) * 

Miok)";:' 

1 k 
+ ( i w) ^ 

-adok);*- 
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(u w) = (u k) , 



+ (i w) " 
1 



, u w 
(i k . , . 
1 k 



This is the two-degree-of-freedom surface 



w u ., , u w 

(uw) = (uk) ^ -r (iw) . -(ik) . ^ 



We shall next consider the slope vector of such a hyper-surface. We take 
partial derivatives with respect to one of the variables, say u, and get 

(u V w 



•)u 


~" 


(u j k . . 


vw 
•'' u j k 

, u' \v 






(i V k . . 


•>i k • ■ • 

U' V 




* 


(i i w . . 


■)i J ••• 



U' V w 
(N-I)(ijk. . .). .^ 



In this, set u = 0. 



(Ovw. . .)^ =(Ojk. . .)J.l 

. 0' w 

M-k. . .). ,^ 



, 0' V 

Mi j w . . . ) . ^ . . 



^ 0' V w 
(X-l)(ijk. . .). . ^ 
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(Ovw. . .)„ = (Ojk. . .)„jk 

0' 
since aU other terms vanish by virtue of =0. 

This result is analogous to the one obtained for boundaries of ordinary 
surfaces; it says that the slope anywhere on a boundary is a function only of the 
slopes at the "ends" of the boundary, and are otherwise independent of the 
bovmdary shapes. Slope continuity across botmdaries is a consequence. 

The hyper-surface equation just developed is defined by ordinary curves, 
or single-degree-of freedom boundaries; we can also write a hyper-surface 
equation for N degrees of freedom, defined by boundaries with N-1 degrees of 
freedom. We shall exhibit the result for N = 3: 

u 

(U V w) = (i V w) . 

+ (U j w) J 
w 

+ (U V k) j^ 
U V 

u w 
-(ivk).^ 

, . , , V w 

- (u 3 k) . ^ 

' 1 J k . 

The proof that this space contains, for example, the boundary subspace 
(0 V w) follows the preceding proofs in principle and will not be carried out. 
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SECTION V 
SURFACE NORMAL VECTORS 

5. 1 GENERAL SLOPE CONTmUITY CRITERIA 

The surface normal vector furnishes a convenient mechanism for the in- 
vestigation of general criteria for continuity of surface slope across boundaries 
between surface patches, tt will be seen that the continuity conditions already 
established are much stronger than are necessary, but that they are expedient. 



Put 



U = fx y z ] 
L u u uJ 

W = [x y z 
'•www' 



for the tangent vectors of a surface patch at some point. Let us assimie that 
another adjacent surface has a common boundary curve along u = constant, w 
varying, so that W is common to both patches. Let the parameter for this second 
patch be v, and for its tangent vector put 



V = X y z I . 
L V V vJ 




COMMON 

BOUNDARY 

CURVE 



/: 
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The two surfaces will be continuous in slope across the boundary at the point in 
question in case the three vectors U, W, and V are coplanar there. 

The surface normal vector for one surface is 

N = U X W. 

If V Is pezpendlcular to N, 

V- N=0. 

But then in this case, the three vectors U, W and V are coplanar since they are 
all perpendicular to N. 



In detail, this gives 

N = UxW= fjjj] 

L X V Z J 



X y z- 



a vector whose components are the familiar Jacobians, and 



N 



'VN 



fx y z I 
L V v vj 



X y z 

V V V 



X y z 
u u u 



X y z 
w w w 



(This is the so-called "scalar triple product" of the vectors.) Thus the vanishing 
of the determinant of the matrix of the three tangent vectors is the general condi- 
tion for slope continuity between two patches, at any point on their common 
boundary. 

This equation also shows that we may have slope continuity of surfaces 
even though the curvilinear coordinates of the two surfaces are not slope - 
continuous across the boundary'. 

If the tangent vectors U and V are equal everywhere along the boundary 
curve, the determinant is sure to vanish; similarly if the tangent vectors U and 
V are scalar multiples of one another, even when the scalar multiplier is a 
variable quantity. 
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GENERAL CONSTRUCTION - TANGENT SURFACES 




Suppose that a siirface A already exists, defined by the parametric vector 
equation 



A = [x{4> e) y{<P B) ^(<Pd )■] 



Let a be a curve on the surface; it is always possible to write the vector equation 
for a in either of two forms: 



a((/. )= [x( <^ ) y( <^ ) z( <^ )] 



a( 0)= [x( 6) y(0 ) z( )] 



Suppose we wish to attach a surface B to surface A, in such a way as to make 
curve a common to both surfaces, and suppose furthermore that we wish to 
maintain slope continuity across this mutual boundary. 

We shall consider curve a to be the boundary (Ow) of the B surface. We 
are at liberty to design, arbitrarily, a projection of the other three boundaries, 
(uO), (ul), and (Iw). Say for example that we design these curves in the xy pro- 
jection. Then the curves represent the x and y components of their coordinate 
vectors. 
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We are now ready to obtain the missing z component of the tangent vector 
across (Ow). We first ccmpute the surfiace normal to A along curve a. For this 
purpose we can use any one of the expressions 

{i>9)^ X (<l> 9)g 

or (4,6)^ X a( e )g 

or a(«^ )^ X {4>6)g 

Each e:q)re8sion yields a surface normal N; the three results are identical. 
We can evaluate this surface normal vector at any point on (Ow) since we have a 
correspondence between w and the variables ^ and d ■ 

We also have the equation 

(Ow) N^ = 0. This is the familiar condition for surface slope 
continuity. 

Let N = [a b c] after ev^duation at w. 

Then the equation becomes 

[x(Ow)^ y(Ow)^ z(Ow)^j 

b 
c 

We alrea^ have the x and y components of this equation, and can solve for the 
z component: 

z(ow)^ = ^ ^<°^)u ^ ^ y<°^>u . 

-c 

This z component has a magnitude that ensures that the complete vector 
(Ow) is coplanar with surface A at w. Hence (Ow)^ is the desired tangent vector 
of surface B across (Ow). 
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We next find iba z tangent vector components z(00) and z(01)^ from (Ow)^, 
and use them in the equation 

1 



z(Ow) 



z(00) F„w + z(01) F.w. 
^ 'u u 1 



Here the 1 superscript indicates that this is an intermediate result; it is the 
intrinsic boundary tangent vector for the F-type surface, and does not yet match 
the z(Ow) vector ftinctiwi obtained from the A surface. 

Accordingly, we must add to the F-type surface a G-type tangent vector 
correction surface, so as to make the combination have the desired slope along 
(Ow). 

This G-type correction surface is, as we have already shown, 



(uw)" = (iGqU G^m] 







Ow 



uO^ 



w 



Iw 



-00 

1 

-10 



uw 

c 
uw 



ul 



-01 



-11 



V 




1 


c 
uw 




V 


c 

uw 




G w 



The superscript c indicates that this is a correction surface. 

Slope correction is necessary only along the boundary (Ow); we can enter 

the value for Ow*' in the matrix, but the other entries must be looked at in detail, 
u 

We have, for the slope correction across (Ow), 

1 



Ow 



Ow 



Ow 



These latter two quantities have already been found for the z component, and so 

c 
Ow is known, 
u 

Consider lw° . This is at a free boundary, (Iw), remote from (Ow), and 

c 

we can set it equal to 0, Then 10 and 11 are both zero also. On the other 

hand u0° and ul° are connected to (Ow) at (00) and (01) and so we must specify 
w w 

them in such a way as to satisfy the conditions at these points. Elsewhere, they 
too are arbitrary. 
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We write the G function expression: 



(uw) = [l GjjU G^u] 



uO 



Ow -00 
u uw 



ul 



-01 



uw 





1 ' 




GjjW 




G w 





(we omit the c super-script temporarily. ) 

Performing the first multiplication, we have 

(uw) = fow G.u I uO -00 G.u I ul -01 Gu] 
* ' I u I w uw > w uw J 



GjjW 



G w 



Consider the element uO - 00 G u In the row vector. Since uO is 
w uw o w 

arbitrary, it would be convenient to choose it so as to make the entire vector 
element vanish. We therefore write 

uO = 00 G„u. 
w uw 



Then 



00 = as it should, and 
w 



since uO = 00 G.u, 
wu uw 



00 = 00 as it should, 
wu uw 



Similarly, we may set 
ul 



01 G„u. 
uw 



The result of these choices of uO and ul is to reduce the G equation to 

w w 



c c 

uw = Ow G.u. 
u 



This represents the correction surface z component which must be added to the 
z component of the F surface in order to obtain slope continuity along (Ow) between 
the given A surface and the designed B surface. 



/ 
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5. 2 ADJACENT- PATCH SLOPE CONTINUITY 

ft is sometiines desirable to define the boundary curves for two adjacent 
patches so that at the junction between the curves the tangent vectors have the 
same direction but are of different magnitudes. This Is particularly useful when 
the boundary curves are parametric cubics, because then the magnitudes of the 
tangent-vectors at the end points control the behavior of the curve segment. 



^ 01 




As a specific case, consider the boundary (Ow) common to two patches^ let 
the tangent vectors at (00) and (01) for the first patch be 00^ and 01 ^ , and let 
the tangent vectors for the next patch be 00 and 01 

If the tangent vectors have the same direction, they are scalar multiples of one 
another, 

or 00 = m 00 , 

u <p 

and 01 = n 01 , . 

u (p 

Suppose that the patch (<^ w) already exists. We need to obtain an appropriate 
e^qiression for (uw) so as to match surface slopes across (Ow). By the results 
of the preceeding article, we can accomplish this in very general ways, but in 
our present case let us make a special requirement on the tangent vectors: let 
us assume that everywhere across (Ow) the tangent vectors have the same 
direction, and differ only in magnitude. Then for any w, 

Ow = ^ Ow^ 
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where X is a scalar. We know that X takes the value m at w=0, and takes the 
value n at w » 1, and we conclude that X therefore must be a scalar function 
of w. 

Witii this relationship between (Ow ) and (Ow^x) the vector cross product is 
always the null vector: 

(Ow)^ X (Ow)^ = [o o] 

Hence the scalar triple product of (Ow) , (Ow) , , and (Ow) vanishes for any 

u *p w 

(Ow) . This ensures that the two surfaces will be continuoJis in slope across 
^ 'w 

(Ow) for anj[ shape of (Ow) and for anj^ X = X (w) that has the proper behavior 
at w = and w = 1. 

We could for example take 

X(w) = m(l-w) + nw. 

This is a linear variation of X with respect to w. It has one disadvantage, 
however, in that it introduces un-wanted cross derivatives or twists at (00) 
and (01). In order to avoid this, we might use 

X (w) = m F w + n F w. 

Then the required slope function across (Ow) for the (uw) patch is 

Ow = (m F„w + n F w) Ow, . 
u 1 <p 

We can check to find tiie cross derivatives introduced by this relationship. 

The cross-deHvative is obtained by differentiating with respect to w, and 

yields 

Ow = (m F' w + n F' w) Ow , 
uw 1 9 

+ (m FjjW + n F^w) Ow^^. 

At (00), 00 = m 00 , 
^ " uw ^w 

and at (01), 01 = n Ow, . 
* ' uw <ffW 
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This shows that the X function does not introduce additional twists at the 
corners of the patch, beyond, of course, those already inherent in the (t^w) 
surface. 

If the (uw) surface already exists, defined perhaps by an F-type equation, 

its intrinsic tangent vector across (Ow) is a known function of (Ow), say Ow , 

where the 1 superscript indicates that it is an intermediate result. Then, 

as before, the correction of slope is 

c 1 

Ow = Ow - Ow , The correction surface 
u u u 

is (uw) 5= Ow G„u, 

u 

When this correction surface is added to the original surface, the combination 

will be continuous in slope with the {^w) surface across (Ow). The u and (^ 

curvilinear coordinates of the two patches will be continuous in slope across 

(Ow), but their tangent vectors will be different in magnitude. 

5.3 APPUCATIONS 

Let dU represent a differential vector, so that 

dU = I dx dy dz I , in which 

dx = r— du + , — dw. 
(Ju Ow 

dy = j-^ du + r^ dw. 
•^ u ow 

dz = J— du + j — dw . 

ou ow 

If dU is tangent to a u-varying, w-fixed curve, these become 

dx = r— du = X du 
du u 

dy = y^ du. 

dz = z du, since dw = 0. 
u 

Hence dU = x y z du. 
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Similarly, if dW is a differential vector tangent to a w-varying, u-fixed 



curve, 



dW 



= X y z dw. 

L WW w J 



The normal differential vector at a point of the surface will be given by 
the vector cross product: 

dn = dU X dW 



du dw. 



The determinants that comprise- the elements of the vector are the 

jabobians J , J and J so that we may write 
X y z 





y z 

'u u 




Z X 

u u 




X y 

u •'u 






y z 

w w 




Z X 

w w 




X y 

w'w 





dn 



= [j J J 1 du dw. 

L X V z J 



X y z • 

The magnitude of dn is equal to the differential area of the elemental parallelo- 
gram described by dU and dW. This magnitude is 



dn I = V dn • dn (or -^J 



M M 



T) 



= du dw VJ 2 + J 2 + J 2 • 
^ ■XT y z'' 

From this, we can construct an algorithm for finding surface areas of patehes; 
we simply perform numerical integration of the expression 







J2+J2''"J2 du dw. 



Again, if N is the unit normal vector to the surface at a point, then 

dn = N I dn , from which 
dn 

hi 



N = 



J J J 
X y z 

s s s 



where 



V J 2 + 



2 -f J 9 + J 2 
x"^ y- z^ 
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The quantities _x, _£, _£ are the coordinate components a, b and c 

s s s 

of the surface normal. 

If the surface Is to be manufactured by milling with a ball-end cutter of 

radius R, then the cutter-center vector Fx y z 1 is related to the surface 
' L c c c J 

vector [x y zl by the simple expression 

k ^c ^1 = ["y ^1 * ^ [*''*=1- 

This describes a "parallel" surface spaced a distance R away from the 
designed surface. 

The normal vector can also be used to calculate volumes enclosed by 
surface patches and planes, as follows. 

Suppose we wish to calculate the volume contained between a surface 
patch and the xz plane. 



dA 



y 

y 
(I 1 — I — I 1 



xz plane 



We can imagine the volume broken up into a number of slender prisms 
whose axes are all parallel to the y axis (and perpendicular to the xz plane. ) 
The area of the base of one of these prisms is the projection of the small 
element of surface area, or 

dA = J du dw. 

y y 
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The x'oUune of tliis prism is 

d\' = \'dA = vJ ciu d\\ 



\- = J J yJ _ du d\\ . 
u w 
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CORNER TWBT VECTORS 



6.1 THE QUASI-SPHERE 

It is possible to choose a parametric cubic that very nearly approximates 
a circle for one quadrant. We shall go into detail about this shortly; intuition 
suggests that similarly we ought to be able to construct an approximation of an 
octant of a sphere by means of a bicubic surface, bounded by these approxi- 
mations to circles. 

For the circle approximation, let us assume that we will be content to 
make the quasi-circle pass through a point on the true circle at u = 1/2. (This 
is not the best possible approximation, but it yields quite good results and the 
arithmetic is simple.) 

We shall assume a circle of unit 

radius, centered at the origin, with 

end-point values of the parameter u 

as shown. The tangent vectors are 

symmetric , but have yet undefined 

magnitudes . 



' 


£ 


1/2 


u = 0^ 


\ 




\ 


u = 1 




' 


-a 



We have 



I u u u 11 



M 



Where the column vector on the 
ri^t represents the end conditions 
for the curve . 
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When u = 1/2, this becomes 
x= ^ [: 2 4 8 



12 4 8 



2 

-3 



1 



1 

-2 

1 



1 " 




" " 


-1 




1 







a 











( - 2 + a + 6 - 4a ^ 4a) 



(4 + a) 



y~2 



a - 8 X - 4, But at u = 1/2, x = ~^ since it is a point on the circle 
(by symmetry, at — ) . 

Hence a = 4 (v 2 - 1). This is the required magnitude of the two tangent 

vectors at u = and u = 1. (Calculation reveals that the quasi-circle has a 

radius of about 1. 00016 at u == 1/3 (at — or 30') so it is a good approximation.) 

6 

We now establish a coordinate system for the sphere, and show its 

boundary curves . 




uO = degenerate curve 
boundary 
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The boundary curves Ow, Iw, and ul are all unit circles; the boundary uO 
however is a degenerate circle, and appears as a point. 

We shall first investigate the z component of the uw surface vector 



z(uw) = F u F u G u G u 
^ ' L 1 1 . 



00 


01 


00 

w 


01 

w 


10 


11 


10 

w 


11 

w 


00 
u 


01 

u 


00 
uw 


01 
uw 


10 
u 


11 
u 


10 
uw 


11 

uw 



F^w 



F w 



G^w 



G w' 



Tf u F u G u G ul 
L 1 1 J 



1 








-a 


1 








-a 



























F w 




F w 



G^w 



G w 



When we perform the first multiplication, we obtain 



5{uw) = [f u * F u I j I -a(F^u -^ F^u) 



F^w 



F w 



GpW 



G w 
1 



But F u - F u = 1, by virtue of the definition of the F functions. 
Hence 



z(uw) =1 



F w 




F w 



G^w 



G w 



F w - a G w. 
1 
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The curvilinear coordinates for w constant thus yield constant z; this implies 
that z is independent of the other variable u, and the w curves are plane curves. 
They must of course be quasi-circles . 

We have obtained the value of the number a by investigating a unit circle. 
For a circle of radius r, the tangent vector magnitudes must be proportional txi 
r , or equal to ra. We can find these radii for various values of w from either 
the quasi-circle y(Ow) or .\(lw). 



We have 
x(lw) 



10 11 10 



11 



-F^w^ 


F w 
1 


G^w 


?l'\ 



and 



y(Ow) 



00 01 00 



01 



F w 




F^w 



G^w 



G w 



in either case, 



r(0w) 



1 a 



'\''' 


F w 
1 


^"^0^^- 


.S"-. 
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For w fixed, the x and y coordinates of a quasi-circle are given by 



F^u F^u GgU G^u 




r 
ar 

J 



and 



V = rFuFuGuGul 
■^ L 1 1 J 



r 




-ar 



where r Is a function of w, shown above. 



But 



and 








r 






= 


ar 




_ 





r 














-ar_ 







1 a 

2 

a a 



1 a 





o 

-a -a' Oj 



Ff^w 



F^w 



'^o^^- 



L^i-j 



^)^^- 



F^w 



G^w 



G^w 



We obtain these last results by simply writing rows in the 4x4 matrix that 
correspond to elements in the vector of the left side. \Mien we comliine 
results , we have 
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ic(uw) = Ff u F.u G u G uj 



"0 








■ 




rvi 





1 


a 












a 


2 
a 
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y(uw) = [FjjU FjU GjjU G^u] 



"o 


1 


a 


' 




[v] 
















F w 

















-a 


2 
-a 







.V. 



The equation for z(uw) has already been shown. 

The striking thing about the B matrix as it appears in these equations is 
that is has non-zero entries in the bottom rl^t partition, 



that 
and 



00 



uw 



10 



01 



11 



x(00^^ = "^ 

y(io^^ = -a' 



By comparison, we see 



These are the croes-derivatives at the corners 00 and 10. All other cross 
derivatives vanish. We shall refer to these cross -derivatives as "twists" of 
the surface; uw is the twist vector at a generalized point on the surface. 

6.2 THE EXACT SPHERE 

A bi-cubic surface cannot fit a si^ere exactly, and it would be interesting 
to see whether by an appropriate choice of F and G functions other than cubics , 
an exact equation can be constructed. 
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The functions 



and 



^" 



2 TT 

F u = sin — u 



resemble the cubic functions in shape; they satisfy the conditions 



F - F 1 - 




F 
1 



F 1 = 



as well as 



^0« 



F^O 



F^l ^ 



^1^ 



1. 



We shall choose these functions, and determine appropriate G functions 
so that the equation 

F^u F^u G^u G^u] UO) j 

(1) 

represents an exact circle, and not the approximation of the last section. 
We already have the well-known parametric equations for a circle; 



X ■" sin au 
y --- cos au 



where a - — 
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u = 



y 


a 


"\ 


" - \ 





We can compare the x equation with 
x(u) 



cos" au sin au G u '^i'-' 



= sin au - a G u 

2 
Then sin au =- sin" au + a G u, whence 



S" 



(sin au - sin au) . 



Again, 



r 2 I 2 j 1 

y ^ cosau cos au 1 sin au G u G u j 



1 





-a J 



whence 



cos au -- cos au - a G u 



G,u -- - (cos au - cosau). 
la 
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We can easily verity that these G functions satisfy the same stipulations 

as the cubic G functions: 

G = G 1 - 


GO = G 1 =0 



S° - ' S^ 



G - G 1 



1 . 



If we now use the same boundar>' value matrices as were used in the 
pre\'ious case of the quasi-sphero , but witli a -= -^ tlirougliout . we obtain for 
the z component of the surface vector 
7. (u\v) F w - a G w . 

This is. witli the new F and G functions 

2 2 

/.(u\v) - cos aw - (cos" aw - cos aw) 

= cos aw. 



.\s before, this shows tliat th^^ z coordinate of the surface is independent 
of u; the w curves are plane cur\-es . and they are indeed circles. Tlieir radii 
are "av-en bv 



.x(lw) 



10 11 10 11 

w w 



^o''^ 



F^w 



S"- 



c; . w 



F w1 


F,w 



G w 





/^l^^J 



F^w - aG^w 

2 2 

sin aw - (sin aw - sin aw) 

sin aw. 
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Then 



x(uw) = If u F u G u G u 





r 

ar 

L J 



and 



y(uw) 



r(F u + aG u) = sin aw sin au. 



F^u F^u G^u GjU 



r 





-ar 



= r(F^u - aGjU) 

2 2 

= sin aw (cos" au - cos" au - cos au) 

- sin aw cos au. 

The resulting parametric equations, when collected, are 
X ^ sin aw sin au 
y = sin aw cos au 
z - cos aw 
and these are well-known. 

This demonstrates that the sphere is a special case of the general surface 
equation, provided the blending functions are suitably chosen. 

The F and G functions are by inspection, seen to be linear combinations of 

the linearly independent functions of u, 

2.2 .1 

cos au sm au cos au sm au , 



and this may be taken as an appropriate basis vector. Then 



F^u F^u G^u G^uJ 



2 2 

cos au sin au cosau sinau 



[Mj 
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where the M matrix is, in this case, 



M 









1 


1 








a 





1 


1 
a 














1 
a 








1_ 
a 






Incidentally, its inverse is 



M 



1 





1 








1 





1 











a 








-a 
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RULED SURFACES 



The locus of straight lines connecting corresponding points on two cuirves 
(uO) and (ul) is a surface. The lines are called "rulings" of the surface. 




The equation for a generalized line of the surface is also the equation for 
the surface: 

(uw) = (ul)w - (uO)w + (uO). 
This is equivalent to 

(uw) = (ul)w + (uO) (1-w). 
From this equation we obtain the derivatives 

(uw)^ = (ul)^^w + (uO)^ (1-w) 

(uw)w = (ul) - (uO) 

(uw) = (ul) - (uO) . 
uw u u 

7.1 DEVELOPABLE SURFACES 

A special case of such ruled surfaces is of importance and interest. If 
the ruled surface has the property of being tangent, along the rulings, to a 
moving plane which rolls around the surface, then it may be deformed by simple 
bending and flattened out into a plane. Such a process is called "development" 
of the surface. We call such surfaces "developables" or "wrapped surfaces". 
(A sheet of paper can be wrapped around the two curves (uO) and (ul) to form 
the surface. These surfaces are also known as "convolutes".) 
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The tangent-plane condition can be established by showing that the tangent 
vector at a point on (uO) , the tangent vector at a corresponding point on (ul) , 
and the tangent vectors along the line joining these points, are all coplanar. We 
need to form the scalar triple proAict of these vectors , and show that it 
vanishes. 

The tangent vectors in question are (uO) , (ul) , (uO) and ("^)^- 

First observe that, for a ruled surface, 

(uw) = (ul) - (uO). This tangent vector is independent of w, so 
that (uO) = (ul) . Moreover, the vector is simply the line segment joining 

WW 

the two points, as ml|^t be expected. We can write, for the scalar triple 
product. 



(UO)^ 

(ul)^ 

(ul) - (uO) 



= 



where the notation represents the determinant of the matrix of the three (row) 
vectors. If the determinant vanishes for all values of u, the surface is 
(tevelopable . 

The preceding describes an analytical test to ensure that a ruled surface 
is developable. We shall now describe a construction that wUl enable us to 
define a ruled surface by means of two space curves. Suppose that the two 
space curves are defined by vector functions of two different parameters , 
u and <f) 
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The scalar triple product is 

(uO) 
u 

(<Pl) - (uO) 

If we consider u the independent variable, the equation enables us to find 
4> for any value of u; this value of ^ determines the point on (c^l) which corre- 
sponds to a point on (uO) , so that the line joining these points is coplanar with 
the tangent to (^1) and the tangent to (uO) . We allow u to vary, and obtain 
related <i> values; these values of <^ enable us to evaluate the components of the 
vector (<^1). These components are the same as the components of the desired 
(ul) vector. 

Provided we remember that the (unspecifiec^ functicms of u and </j are 
different, the symbolism 

(ul) = (<^1), <^ = <^ (u) 
represents the statement that the Fx y z ] vector Is the same for both. 

With the correspondence established between points on the two curves, 
we can write the equation for the developable surface , 

(uw) = (ul)w + (uO) (1 - w) . 
This is the ruled siirface equation, but with a special relationship between 
curves (uO) and (ul) . 

7.2 PLANE/SURFACE INTERSECTIONS 

The general sur&ce equation can be cast in the form 

T 
uw = UBW where U and W are vector functions of u and w 

respectively, and where B is a square matrix describing the boundary curves. 

For example, we might be dealing with the first F-type surface equation. 



■-[• 



1 F u F 



■"] 




Ow 

Iw 



uO 
00 
10 



ul ' 




-1 


01 




F^w 


11 




F w 



in which these vectors and the matrix are explicit: 
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Since uw is in reality a vector consisting of an x, a y and a z component, 
there are three B matrices \vhich we can call B B B^. 

We wish to find the intersection of this surface with the plane 

ax + by + cz + d = 0. 

We can substitute x = UB W^, y = UB W^ and z = UB W into this equation, 
X y z 

and write the result in the form 

U faB + bB + cB 1 W + d = 0. 
'■ X y z J 

It is permissible to interchange the order of multiplication from aU, bU, 

cU to Ua, Ub, Uc because a b and c are scalars. In this form, the sum 

faB + bB + cB 1 = S, a square matrix function of u and w, and 

i X y zJ 

USW^ = -d is an equation in the two variables u and w. If w is 
assigned a fixed value, there results an equation in u which when solved will 
yield a point on the Intersection curve of the surface with the plane, (provided 
of course such a point exists for the chosen value of w.) 

If the surface in question is a bi -cubic , the matrix S is no longer a 
function of the variables u and w, but consists of constant elements. In this 
case the above procedure reduces to the solution of a series of cubic equations 
in u, where the coefficients of the cubics are determined by successive fixed 
values of w. 

In any case. If the spacing of the w values is close, the old value of u 
just previously determined for a particular choice of w can appropriately be 
used as a first trial solution for the new value of w. Algorithms for the 
improvement of this initial trial value of u are not difficult to construct, and 
will not be discussed in detail. 

If the plane is given by, say, the equation 
x + d = 
the solution procedure is unaltered. Not much simplificatton results. 
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RATIONAL POLYNOMIAL FUNCTIONS 

8. 1 BOUNDARY CURVES AND BLENDING FUNCTIONS 

Two kinds of curves have for many years been traditionally used in air- 
plane lines design — cubic polynomials, and conies. Unfortunately, each of 
these curve forms for itself has certain drawbacks. In the parametric form, 
for ordinary cublcs, the entire shape of a curve segment is governed by end 
tangent vectors. Sometimes these end tangent vectors lead to unwanted hooks 
and bulges in the curve segments. On the other hand, conies, although more 
benignly behaved, cannot by their very nature yield curves with points of 
inflection. Yet such curves very often exist in aircraft shapes — as for in- 
stance in the case of wing fillets. 

Because of these short-comings, a new curve type has been developed. 
It is based upon rational polynomial functions. It contains both conies and 
ordinary cubics as special cases, and provides a great degree of generality 
and flexibility. 

We start by establishing the form of this function. 

Let v be a vector, so that for example v = [x y z 1] or v = [x y 1] or 
V = [x 1] . The first of these can be thought of as the vector (or matrix) of 
coordinates on a space curve; the second is the vector of coordinates for a 
plane curve, and the last is the vector of a single coordinate. Since this last 
vector yields the most general case, we shall begin with it, and show how one 
might evaluate a set of numbers in a matrix to define each of the parametric 
coordinates of a curve. 

The product of v and a variable scalar w is wv = [wx w] . Here both 
wx and w are cubic functions of a parameter, v, and obviously 

wx 
X = . 

w 
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This 1b the ratio of two cubic polynomials (hence the term "rational func- 
tion"). We can represent the two cubic polynomials by the matrlc equation 

wv = [u3 u2 u 1] A. 

Since wv •: [wx w] , the matrix A must consist of four rows and two columns 
of constant coefficients. We now proceed to show how these numbers may be 
found so as to define a coordinate of a curve. 

We shall be interested in the end-point coordinates of the curve at v :> 
and V •= 1. These coonUn^es are v^ = [x,, 1] and Vj^ = [x 1] respect- 
ively. A tangent vector aiqrwhere on the curve is clearly 

V' = [X' 0] 

where the prime mark means differentiation with respect to the parameter v. 
The tangent vectors at v » and v « 1 are therefore v ' = [x ' 0] and 
V ' = [X ' 0] respectively. 

Now differentiate both sides of 

wv ' [tt3 u^ u 1] A and obtain 

(wv)' - [3u2 2ul 0] A. 
Substitution of u ' and u = 1 Into these two expressions yields 

'o r 

1111 
10 
3 2 10 



*0^0 



*1^1 



<^oV 



(Wj Vj)' 



A. 
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The 4x4 matrix on the right has an inverse, and we may write 



1 


-1 


'"o^o" 


1111 




W^Vj 


10 




(*oV 


3 2 10 




.<*'l ''l)'. 



2-211 
3 3-2-1 
10 
10 



^0 


^0 






"l 


'l 






w ' 
*0 


^0 


*^0 


-o' 


w ' 

1 


^1 


^"l 


^l' 



The square matrix inverse Is constant and always the same, and reawwars in 
the algebra so often that we shall henceforth call it the matrix M. 

The matrix equation can be factored and rewritten in the form 



A = M 



^0 













^0 





^1 










\ 


-o" 





*0 







\ 





w; 





*1. 




."l 



The right hand matrix of v's represents the desired end conditions on the 
curve. In our present case, it is of course a 4 x 2 matrix. 

The middle matrix is 4 x 4 and contains the four numbers [w^ w^ w^' vi^']. 
Any arbitrary set of four numbers inserted into this matrix will serve to define 

a unique pair of cubic functions of the parameter u, from which x can be found, 

wx 
by using, as we have said, the ratio x = ~- . 

Instead of picking these four numbers arbitrarily, however, we shall 
Impose further conditions on the curve until enough conditions are imposed 
to define [w. w w ' w '] uniquely. 
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We begin by Introducing desired second derivative vectors at the end- 
points; these vectors are clearly 

Vj," = [Xjj" 0] and 



Vj" = [Xj" 0]. 

(Incidentally, in the case of vectors v' = [x' y' 0] , v" = [x" y" 0] if 
the determinant of the matrix fx' y'l vanishes, the curve will have a 

[x" y'J 
point of inflection at v. If the determinant Is positive, the center of 
curvature will lie on the left as one proceeds along the curve; if the 
determinant Is negative, the center of curvature lies on the right. If 
two curve segments have equal v' and v" at a junction, they are contin- 
uous both In slope and curvature at such a junction. ) 

When we take second derivatives of both sides of 
wv = [u3 u2 u 1] A, we obtain 
(wv)" = f6u 2 0] A. 
At u = 0, this is 

(^oV" = [0 2 0] M [wqVjj 

w^v^ 

|_W ' V + W V ' 

But (w„ v„)" = w„" v„ + 2 w' v' + w„ v„" and solving for w„ v„", 
*0 0' 00 00 00 00 

*o \" = <^o V" ■ *o" ''o ■ 2 ^o' ^o'- 



Now (w V )" = [-6 6 -4 -2] 



^0^0 



w^V^ 



*o' ''o -^ "o ^o' 



'-*! \ * ^ ^IJ 



/^ 
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Furthermore, since in general 
(wv)'' = [(wx)" w"] , 

the quantity w" is the second component of the vector of (wv)" and therefore is 
associated with the last column of the matrix 



^^'o ^0 






"'l ^1 






^o'^o 


+ 


^^0 ^o' 


w; v^ 


+ 


-1 ^: 



But the last components of v and v are both 1, and the last components of v^' 



and v ' are both zero. Hence 



[-6 6 -4 -21 



IV 



1 



We now can write, b}' combining results 
[-6 6 -4 -21 



^^'O ^'O 



^^0 <^'o - \^ 

Wq' (Vj^ - v^) + Wy V 
[-V (V, - v^) 4 w^ v^' 







-2 w^. v^' 



6 w 



'•i(^'i-V--'%V--"i'^-V-^"i^'i 



IJ 



Collecting 



■0-0 •■o(-^V)-^^^i(«^-V--^i') 



w' v„" = w. 



;(-2v^') . w^^' (-2 (v^ - v^)). 
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We now restore this last ejqtresslon to matrix form: 



Vo" - f*o "i *o' "I'J 



-4v. 



6 (Vj - Vq) - 2 Vj' 



-2 V. 



-2 (V, - V J 



In our present case, the matrix on the right consists of a column of numbers 
and a column of zeros. Hence the column of zeros can be discarded, and the 
result Is a 4 X 1 matrix. On tke left, w v " la a scalar. 

Similarly, we can find by analogous algebraic procedures that 



Vi" = f*o *i *o' *i'l 



'6(v„-vj + 2v. 



4vi' 



-2v ' 



Then, writing a matrix equation, we have, so far, 

where [ P | Q] represents a 4 x 2 matrix consisting of the separate 4x1 matrices 
for w- V " and w v ", written side by side as coliunns. 

We now introduce another condition. Let it be required that the curve 
pass through the point v = [x 1] when u = — 



(This value of u Is of 



course arbitrary. ) 

This condition leads to 



V = -=-[124 8] M 
c 8 ' ' 



*0^0 






W^Vj 






V^O 


+ 


*oV 


"I'^l 


+ 


^v. 



/; 
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By algebraic manipulations similar to the preceding, we can rearrange the 
equation to read 



8v^ = ["0*1%'"!'^ 







-V, 



This is an equation in which v has two components, x and 1. The right hand 
matrix Is a 4 x 2 . Call it the R matrix. Then we can adjoin these matrices, 
to obtain 

[Wq Vjj" w^ v^" 8v^] = [Wg w^ Wp' w^'] [ PQ RJ. 

Now [ P Q R ] represents a 4 x 4 matrix; P and Q are each 4x1 matrices, but 
R is a 4 X 2 matrix. We next transfer w^ v^^" and w^^ v^^" to the right hand 
side, obtaining 



[0 8vj = [Wjj w^ Wg' w^'] /[PQR] 



^o" 











-l" 
























°1\ 






oj / 



The right hand matrix is now a 4 x 4. Provided its determinant does not vanish, 
it has an inverse, and 

[ Wq Wj^ Wjj- w^'] = [0 |0 |8v^] S 

where S is the 4x4 inverse of the matrix. 

Now that [ w w w ' and w ' ] have been evaluated, the curve is com- 
pletely defined, since the rational function — — is completely defined. 



/r 
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8. 2 PLANE CURVES 

With some loss of generality and flexibility, we can have v = [ x y 1] , 
a plane curve. We shall next show that for an appropriate choice of 
[ w w w ' and w. '] , the curve reduces to a conic. 

We have the equation 

wv = [ wx wy w] 

= [u3 u^ u 1] A. 

In this case, A is a 4 x 3 matrix. Now if the top row of this matrix is 

[0 ] , the equation reduces to 

2 
wv = [u u 1 ] A 

when the top row of A has been omitted. A is now a 3 x 3 matrix, and it is 
possible to show that this equation is a parametric form for the general conic, 
expressed as a quadratic rational ftmction. 

For the top row of A, we have the vector equation 



[2-211] 



*0^0 






w^v^ 






*o'^o 


+ 


*oV 


V^l 


+ 


Wj v^'_ 



=[0 0] 



Expanding: 

2 Wg Vp - 2 w^ v^ + Wq' v^ + Wq v^' + w^' v^ + w^ v^' = [0 ] 
Collecting: 

Wg(2vQ+VQ') + w^(-2v^ + v^') + w^|' (Vj|) + w^'(Vj^) =[0 0] 
In matrix form. 



[0 0] = [Wjj w^ Wp' w^'] 



2 v„ + v„' 




-2 v + V ' 
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The matrix is a 4 x 3. 

Now we can adjoin a column to the matrix and an element to the vector, 
and write 



f Wp] = [Wq w^ Wq' w^'] 



then 



Wp Wj Wq' w^' 



] - [0 Wq] 



2^'O^V 



-2 v^ + v^ 



- ^'o ^ ^'o' 


1 







^0 





^'l 






-1 



If the indicated inverse exists, then a solution can be obtained in terms of w^. 
Furthermore, w can be set equal to 1 arbitrarily. 

The matrix has an inverse in case the determinant 

-2 V + V ' / 0. 







1 



As a test, construct a conic with end conditions 

Tool" 
I 1 1 1 
; 1 







L^l 



i 1 
L 
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2^0+V 



-2 V + V ' 



[ 2] + [ 1 0] 
[ -2-2-2] +[100] 
1 
111 

We first test to see whether the determinant vanishes: 
-1 -2 -2 



10 2 

-1 -2 -2 

1 

111 



1 
111 
will have an inverse 
The matrix is 



-1 



•1. Hence the augmented matrix 



1 


2 


1 


1 -2 


-2 








1 





1 1 


1 


0_ 



and its inverse is 






1 





2 





-1 


-1 


-1 








1 





1 


-1 


-2 


-2_ 



then 



[ 1 w^ Wp' w^' : 



[ i; 



10 2 
-1 -1 -1 

10 

1 -1 -2 -2 



whence [ 1 w w ' w ' ] = [1 -1 -2 -2 
The conic equation is 



wv = [ u3 u"^ u 1 ] 



2-211 

-3 3 -2 -1 

10 

10 



1 







-1 







2 


1 




-2 


-1_ 








1 


1 1 


1 


1 





1 


0_ 
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[ u^ u2 u 1 ] 



1 -ll To 1 
1-1-2 1 111 

-2010 100 

1 oj I 1 0_ 



= [ u^ u2 u 1 ] 









o' 


-2 


-1 





1 





-2 








1 



= [u^ul] 



-2 -1 
10-2 
1 



[wx wy w] = [-2u +u |-u J -2u + ij 



y = 



-u 



1 - 2u 
asymptote at 1 - 2x = 0, x 



1 
2 



2x 



The curve is hyperbolic, with an 



It is always possible in aU of the foregoing to set w^j = 1. This is 
because aU equations are homogeneous. It Is never possible for w^ = 0, 
since this leads to certain degenerate cases. 

We remark in passing that when 

[W(j Wj Wjj' w^'] = [1 1 0] 
the equation reduces to the ordinary parametric cubic, given by 
V = [ u3 u2 u 1 ) M 



and w is constant and equal to 1. Hence the rational polynomial fimctions 
contain as special cases all conies, ordinary cubics, and of course therefore 
straight lines and circles. 
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Their use as boundary curves for surface patches is obvious. They 
maintain tangent vector continuity between adjacent patches; indeed, if the 
F. and F^ functions are constructed as rational functions, we can establish 



the F, fu 


nctiot 


I 






Fj{t 


, "3 


') = , 


Su" - 3u + 1 


This function has the end conditions 




'^0 ■ 




"o 1 " 






\ 




111 






\' 


- 


10 






\' 




10 






v 











.-l". 










Since Vjj" = \v ' = *( 1 0] and v^" = \v ' = [1 0], the 
curve has a point of inflection at u = and u = 1. Hence its use insures 
curvature continuity across boundaries between patches, po-ovided of course 
the boundary curves have similar curvature continuity at patch corners. 

The cure is symmetric. Furthermore, we can put 



Fo(u) 



1 - F^(u) 



and obtain directly the F. ftinctlon, another cubic rational function, with 
similar properties to F . 

8.3 AN EXAMPLE 

We shall work out the equation for the Fj(u) blending function with the 
customary stipulations that 

F^(l) = 1, F^(0) = F^'(O) = F^'(l) = 0, 



and with the two additional stipulations that F "(0) = F "(1) = as well. 
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This blending function will give both slope and curvature continuity across the 
common boundary between two contiguous patches. The end-conditions are, 



for V = 


[Fo(u) 1], 
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The matrLx [ P Q R 1 


- 




^^'o' 
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_0 
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-6 





4^ 








(i 


4 


4 








-2 





1 








1 


-2 


-1 


_j 


L ', 







obtained bv direct substitution in the given form. 



Its inverse is 





-1 





-4 




1 





-4 







_2 


2 


6 


6 




2 





-6 
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Now set v„ 



[^.] 



for symmetry. Then, since u = 



[Wq w^ Wj^' w^'] = [0 8 v^] S 



[0048] 



0-1 0-4 

10-40 

-2266 

2 0-60 



[8 8 -24 24]. 



Now if w = 1, instead of 8, the equation becomes 



[Wq w^ w^' w^'] = [1 1 -3 3] 



we have 



A = M 



^''o^o 






"l^ 






"o' "o 


+ 


^o\' 


'\' \ 


+ 


^ \' 



and substituting the values of [ w w w ' w ' 



2-211 

-3 3 -2 -1 

10 

10 



1 

1 1 
-3 
3 3 



3 2 



Finally, wv = [ wx w ] = [ u u u 1 ] 





= 


'l 










3 







-3 







1 




1 


o" 









3 









-3 









1 





[ wx w ] 



[ u^ I 3u2 - 3u + 1 ] 
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Hence x = 



F^ function is 



wx 
w 



3u - 3u + 1 



F (u) as required. The other 



Fq(u) = 1 - F^{u) = 1 



3u - 3u + 1 



3 2 

-u + Z\x - 3u + 1 

3U - 3u + 1 



This Is seen to be a rational cubic function also. 
8. 4 PLANE CURVE THROUGH A POINT 

The plane curve vector isv= [xy 1], wv= [wxwyw]. Here 



the polynomial denominators in x 



wx 

w 



andy 



wy 
w 



are both the same. 



As before, the end conditions on the curve are contained in the 



matrix 



We wish to cause the curve to pass through some arbitrary point v 
(commonly called a "shoulder point") and it wiU turn out that we shall also 
be free to choose some arbitrary slope at this point. It is important to 
distinguish betwera the term "slope" and "tangent vector". The slope of a 
curve implies that the direction of the tangent vector is known, but the magnitude 
of the vector is not under our control. 

We begin by assuming some value of the parameter u to correspond to 
V . For purposes of illustration, let u « —r- at this point. 
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3 2 
Then, from wv = [ u u u 1 ] A, we write 



-—[12 4 8] M 



*0^0 






*1^ 






v% 


+ 


*0% 


"I'^l 


+ 


^^1 



Observe that we have arbitrarily set w^ = 1. This is harmless, since 
the equation is, as we have observed, homogeneous. 

By multipljrtng^the matrices, combining, and collecting terms as we 
have done before, we achieve the result 



8v^ = fWo \ "o' "l'^ 



4 V - v ' 
1 1 



The matrix on the right is a 4 x 3 matrix; to make it square, so that it 
can have an inverse, we need an additional column. This column can be pro- 
vided by a scaler equation, and the slope relationship will fUrnish this equa- 
tion. 

We first find an expression for the tangent vector at v^. Differentiating, 
we obtain as usual. 



(wv)' = (3u''2ul 0] A 
w' = [ 3u^ 2u 1 ] M 



w« 



y. 
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This last equation comes from the equation for (\vv)' by an argument that 
we have used before; that since wv = [ wx wy w ] , (wv)' =^ [ (wx)' (wy)' w' ] . 

Hence w' corresponds to the last column of [w v 



w^v^ 






which is 



Now at V , we have already set u 



We make this substitution, and 



obtain, from 



wv' = (wv)' - w'v 



[-66-11 



™1 (\ - \^ 



' 







™l' (^1 - ^c> ^ '\ \[ 

When we perform the indicated multiplications, and then collect results 
and restore to matrix form, we have 

4v; = [Wg w^ w^' w^'] --6 (v^ - v^,) - Vq'- 

6(v^-v^)-v; • 

^'c - ^0 



V - V 

c 1 
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We are now ready to introduce the slope condition. We could write 



Jl 



; , but thia would yield awkward results when the slope became very 



djr 

dx X' 

great and aj^oached an infinite slope. Instead, we choose two nunbers a 

and b so that ax' <= by'. It is obvious that these numbers can very appropri- 
ately be a = sin 9 , h " cos 9 where d is the slope angle. Then, for v , 



^y„ 



ax 



This is a scalar equation. 

Now x' is the first component of v' , and y' is the second component 
c c c 

of V . Hence these quantities correspond to the first and second columns of 
the matrix in the equation for v' , respectively. We can write this out in 
detail: 



= [Wjj Wj Wjj' w^'J 



-6<Xj,-x^) - Xj,' 
6(Xj-x^)-x; 



X - X- 

c 



X - X, 

c 1 



-«(yo-V-V 






eorj-y^hy^' 






yc-^o 






^c-^i . 






= [Wjj w^ w 


d' V^ 


'P 

q 

r 
s 



where the brackets on the right enclose the resulting 4x1 matrix (or column 
vector). 



The factor 4 (of 4 v ') obviously drops out of the equation. 
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We now adjoin this to the equation for v . Again we can drop the 8 (from 

c 



8 V }, and obtain 



l\ 0] = ["o *1 *o' V 



4 V - V 
1 1 



The vector on the left consists of four components: [ x^ y^ 1 ] ; the 
matrix on the right is a 4 x 4 matrix. If it has an inverse, S, then we can 
solve for [ w w w ' w ' ] by the equation 

[Wjj w^ Wjj' w^'] = [v^ 0] S. 

These values of fte w vector cause the curve to satisfy the desired 
conditions. 

8. 5 SECOND DERIVATIVE VECTORS 

We have already discussed rational functions for 

v = [ X 1 ] (and of course [ y 1 ] and [ z 1 ] . ) 

In particular, we showed that curves based upon these functions can 
usually have arbitrary first and second derivative vectors at the end-points, 
and in addition can be caused to pass through some shoulder point v^, also 
arbitrarily chosen. 

When the vector v = [ x y 1 ] , the complete generality of the resulting 
curves is somewhat curtailed. We shall investigate the conditions under which 
such a plane curve can satisfy end conditions including second derivatives. 

We have already obtained an expression for the vector 

Before, the vectors w- v " and w v " were actually scalars, since they came 
from v" = ( x" 1 . But since v" = [ x" y" ] , they are each 2-component 
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vectors, and their combination makes a 4-component vector. Similarly P aitd 
Q are now each a 4 x 2 matrix, and their combination is a 4 x 4 matrix. 

We carry w. v " and w v " across the etjual sign, and obtain on the left 
the null vector: 

[0 0] = [Wjj w^ w^' w^'] /[p|q] 
















The matrix on the r i^t is 4 x 4. Kow the condition that must hold, in 
order for there to be a solution for [ w w w ' w ' ] , is that this matrix 
must be singular; the determinant of this matrix must vanish. 

This last remark tells us that v." and v " cannot be chosen entirely 
arbitrarily. However, it is always possible to make the determinant of the 
matrix vanish by the adjustment of any one of the four components of v." and 
v ". Thus if one of the four components is the number a, we can expand the 
determiiiant in such a way as to obtain the equation. 

k a + k B 0, from which a can be found. 

Suisse the matrix is, or has been caused to be, singular. Then, if 



[0 0] = [Wji Wj Wjj' Wj'I S, 

we make it non-singular by an appropriate modification. In some cases, this 
might consist in adding 1 to an element in the top row of S. As an illustration, 
we mig^t have 



[0 Wjjl = [Wjj Wj Wjj' w^'] 



S + 



/? 
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K can be seen that this modification is still a valid equation; if the 
modified S matrix now has an inverse, we can immediately obtain it and solve 
for the w vector. 

As an illustration, consider the end conditions 

1 

111 

10 

10 

10 

V " a 

We plan to adjust a in v " until the matrix is singular , 

We require first the matrix [ P | Q ] , which is given by 



-l' 



-4Vo' 


6(v,-v^) 


6 (v^ - Vjj) - 2 v^' 


4v; 


-2v; 


2(Vo-v,) 


-2(Vi-Vo) 


-2v; 



Substitution of the end conditions gives the 4 x 4 matrix: 



-4 

[6 6) - [2 0] 

-2 

-2 -2 



[-6 -6] + [2 0] 

4 

-2 -2 

-2 



-4 0-4-6 

4 6 4 

-2 -2 -2 

-2 -2 -2 



In passing, we note that if our end conditions had been v " = f ] 
and if also v " = [0 ] , the resulting matrix would be singular, because 
the first and third columns of [ P I Q ] are identical. However, this is not 
our laresent case. 



/ 



96 



SECTION vm 



We now subtract the matrix 



\' 


o" 




1 





o' 





^l" 










a 




= 











































from [ P I Q ] and obtain 

-5 -4 -6 

4 6 (4-a) 

-2 -2 -2 

-2 -2 -2 



By a series of reductions accomplished by multiplying rows of the matrix 
and additions (or subtractions) of rows to remove elements, we can obtain the 
determinant 

(a-1) - 1 

1 - 1 

This implies 



This is the value of a that makes the matrix singular. 



The singular matrix is 



-5 -4 -6 

4 6 2 

-2 -2 -2 

-2 -2 -2 

We make it non-singular by adding 1 to the top left element. The matrix 



-4 -4 -6 

4 6 2 

-2 -2 -2 

-2 -2 -2 
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and corresponds to the vector [ w ] on the left of the equation. 
The inverse of this matrix is 

1 



2 1-63 
-1 3-1 
-1 -1 3 -3 
-10 2 



R. 



Fi 



nally, [Wq w_^ w^' w^'] = [Wg 0] R. 



If we arbitrarily set w = 2, then the required solution is just the top 



row of R, or 



[Wp w^ Wq' w^'] = [2 1 -6 3]. 

With these numbers known, the curve equation is completely defined. We 
obtain it bv substition in the canonical form: 



3 2 
[ u u u 1 ] M 



w v 




w^ v^ 



^^o' ^0 ^ ^0 ^o' 
^' ^1 ^ ^1 \' 



u u" u 1 ] M [002 
! 1 1 1 



3 2 , , 
[ u u u 1 ] 



[0 -6] + [2 0; 
[[3 3 3] + [1 0] 



2-211 

-3 3 -2 -1 

10 

10 



2 
111 
2 0-6 
4 3 3 
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3 2 
u u u 1 1 4 1-1 

-5 6 

2 0-6 

2 

The separate equations for the x and y coordinates of the curve can be 

exhibited: 

3 2 

4u - 5u + 2u 



3 2 

-u + 6u - 6u + 2 



-u + 6u - 6u + 2 



In the foregoing, certain matrices have occurred. These matrices are 
significant ones, and can be written as transformations of the common 
matrix 



^o' 



L^ 



as follows: 










For the conic condition matrix, 








''^o^^o' "" 




~ 2 


1 







-^\-\' 




-2 





1 




\ 




1 










\ 




1 
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For the u = — , v shoulder point condition, 
2 c 



4 ^0 ^ ^o' 



4 v^ - v^' 



4 10 
4 0-1 
10 
0-100 



For the P and Q matrices associated with w v " and w v " 













-4 







'^o" 


6 v^ - 6 v^ - 2 v^' 




-6 


6 





-2 




\ 


-2 v^' 










-2 







\' 


-2 v^ + 2 v^ 




2 


_2 










/l'. 


6v^-6v^.2v/ 




6 


-6 


2 







^0 


4v^' 













4 




^ 


2v^-2v^ 




2 


-2 










^o' 















-2 
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